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Abstract: We generalize the computations of ||] to generate long wavelength, asymptoti- 
cally locally AdSs solutions to the Einstein-dilaton system with a slowly varying boundary 
dilaton field and a weakly curved boundary metric. Upon demanding regularity , our solu- 
tions are dual, under the AdS/CFT correspondence, to arbitrary fluid flows in the boundary 
theory formulated on a weakly curved manifold with a prescribed slowly varying coupling 
constant. These solutions turn out to be parameterized by four-velocity and temperature 
fields that are constrained to obey the boundary covariant Navier Stokes equations with a 
dilaton dependent forcing term. We explicitly evaluate the stress tensor and Lagrangian as 
a function of the velocity, temperature, coupling constant and curvature fields, to second 
order in the derivative expansion and demonstrate the Weyl covariance of these expressions. 
We also construct the event horizon of the dual solutions to second order in the derivative 
expansion, and use the area form on this event horizon to construct an entropy current for 
the dual fluid. As a check of our constructions we expand the exactly known solutions for 
rotating black holes in global AdSs in a boundary derivative expansion and find perfect 
agreement with all our results upto second order. We also find other simple solutions of 
the forced fluid mechanics equations and discuss their bulk interpretation. Our results may 
aid in determining a bulk dual to forced flows exhibiting steady state turbulence. 
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1. Introduction 

The gauge gravity correspondence [|| is one of the key ideas to have emerged from devel- 
opments in string theory. It has already led to considerable progress, and promises to 
teach us much more both about quantum gravity and about strongly coupled field theory 
dynamics. 

The best studied class of gauge gravity dualities relate the dynamics of a particular 
conformal field theory on R 3 ' 1 to the physics of a corresponding gravitational theory on an 
AdSs background. We now have an infinite number of proposed dualities of this nature. 
As quantum field theories may be formulated on arbitrary base manifolds it is natural to 
attempt to generalize these AdS/CFT correspondences to obtain dual descriptions of these 
CFTs on an arbitrary Lorentzian base manifold M%\. 

According to the usual rules of the AdS/CFT correspondence (See for a review), 
the classical phase space of the dual bulk description is the set of regular solutions of the 
relevant bulk equations of motion that asymptote at small z (the boundary) to the metric 1 

dz 2 + ds\ i 

ds\ = -g-M (1.1) 

where ds\ 1 is the metric of Ms^See Appendix |C| for a list of notation employed in this 
paper). In this paper we use the methods presented in ||] (see also |], 0) to com- 
pletely characterize a special corner of this phase space that is dual, under the AdS/CFT 
correspondence, to boundary fluid dynamics. 



1 Together with similarly prescribed asymptotic conditions for all other fields. 
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Through most of this paper, we focus on the bulk dual of the four-dimensional CFT 
dynamics at large N and infinitely strong coupling formulated on a weakly curved four 
dimensional spacetime 2 . Away from possible singularities, the the bulk equations that 
determine dual dynamics in this limit may be derived from a two derivative action that 
describes the interaction of gravity with other fields 3 The equations of motion of all 
such systems have a universal subsector, in which the Einstein frame metric is a solution 
of Einstein gravity with a negative cosmological constant and all other fields are simply 
zero 4 . In this paper, we study this universal subsector subject to the boundary conditions 
(0). 

We are specifically interested in regular bulk configurations whose variations in the 
boundary directions are slow compared to the length scale set by the local energy density 
of the solution. Experience with field theory along with AdS/CFT suggests that dynamics 
of such long-wavelength solutions is effectively described by four-dimensional hydrodynam- 
ics. 5 . Starting from the bulk Einstein equations, in this paper we demonstrate that this is 
indeed the case. More specifically, within a boundary derivative expansion, we construct 
an explicit map from the solutions of Navier Stokes equations (with distinguished, gravi- 
tationally determined values of parameters) on the manifold M3 i to the space of regular 
solutions of Einstein's equations with a negative cosmological constant that asymptote to 
dDD at small z. 6 

The method employed in this paper is a direct generalization of the procedure employed 
in ||] (see also |J). Given a velocity field u^(x) and a temperature field T(x) we note 
that the metric 



ds 2 



-Zu^x^dx^dr - r 2 f(b{x> x )r)u IJi u v dx ,1 dx v + r 2 V ^dx^ dx l 



,„, - 5^(x M ) + u^x^u^x^) 

1 

A 



(1.2) 



1 



(where g^ u (x^) is the metric on M^\) has several desirable properties. First, its asymptotic 
form matches ( |1 . 1[ ) . Next note that lines of constant x^ 1 in ( |1.2[ ) are ingoing null geodesies. 
Tubes centered around these geodesies are locally similar to tubes in a uniform black-brane 
geometry but with a temperature and velocity that varies with x^. In fact, this metric 



2 In this introduction, we will first describe our results assuming no variation in the boundary coupling 
constant, but we will later generalize these results to a prescribed weakly varying coupling constant. 

3 For example, when the CFT in question is J\f = 4 Yang Mills theory, the action corresponds to the S 5 
compactification of IIB supergravity. 

4 A class of two-derivative gravitational theories also admit a consistent truncation to a larger Einstein- 
dilaton system which we will employ to describe later the bulk duals of boundary theories with prescribed 
variation of the coupling constant 

In general, the fluid mechanics approximation is self consistent only when the length scale associated 
with curvatures of the boundary metric ds§ i, and the length scale of variation of the boundary coupling 
constant, are also large in units of the same energy density scale. We assume that these conditions are met 
in the following analysis. 

6 See M p|, p|, p], |i^ p|, p], p|, p| p], p|, pi, M, P, H pi pi pi M pL pi, pi K P, Eil 0, Em 




61, 7|"62[ 6^" S4l f° r a us t °f relevant literature on this subject 



5j |Kj, |5^, |5j, |5£| 
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reduces to the exact uniform black-brane solution of Einstein's equations when u^{x) and 
T{x) and g^uix) are constants. Finally, under suitable assumptions regarding the late-time 
behavior of the boundary-metric, boundary velocity, and the boundary temperature fields, 
this metric has a regular future event-horizon. 

We then generalize the procedure of |l]] to demonstrate that ( |1.2[ ) may be thought of 
as the first term in a systematic expansion of Einstein equations in a power series in yz 
where T is the local temperature of the solution and L the local length scale of variation 
of the temperature and velocity and metric functions in the boundary directions. 

We work out the rules of this perturbation theory to arbitrary order in e, and explicitly 
implement these rules to second order. As in [01], this procedure works only when the 
velocity and temperature fields obey a particular dynamical equation of motion. This 
equation of motion turns out simply to be V^T^ = (the indices (p, v) run only over 
the boundary directions ) where the boundary stress tensor T^ u is a given function of the 
velocity and temperature fields, whose form is determined, order by order in the derivative 
expansion, by the perturbative procedure itself. Of course this equation of motion is simply 
the Navier Stokes equation of fluid dynamics, with particular values of fluid parameters 
predicted by gravity. Our explicit second order solution to this perturbation theory yields 
an explicit expression for the fluid dynamical stress tensor to second order in the derivative 
expansion. Most terms in this stress tensor are obtained by the simple covariantization 
of the flat space second order stress tensor reported in jl|]. However in addition we find 
a contribution proportional to a curvature tensor that the analysis of Q was blind to. 
The existence of such a term was already predicted in [pQ ], and our result is in complete 
agreement with their prediction (including coefficients) providing a check both of our results 
as well as those of |50t| . 

We show on general grounds that the metrics constructed by the perturbative proce- 
dure developed in this paper must transform in a specified fashion under boundary Weyl 
transformations. In particular, Weyl covariance demands that the metric dual to a fluid 
configuration be of the form 

ds 2 = -2Wi u^dx^{dr + rA x dx X ) 

+ [r 2 (W 2 g^u + W 3 u„u v ) + r {W^u u + u^W Au ) + W 5fiu ] dx^dx u . ( L3 ) 

Ax = u.Vux —u x 

where Wi, W% and W3 are Weyl-invariant scalars, W^u is a Weyl-invariant transverse vector 
and Wsuv is a Weyl-invariant transverse traceless tensor and Ax is the fluid mechanical 
Weyl connection introduced in fl5"l| ]. This requirement, while logically necessary is not 
algebraically automatic in our perturbative procedure, and so gives a second nontrivial 
check on our results. Our construction passes this test; in fact we are able, in §||, to rewrite 
the bulk metric dual to second order fluid dynamics entirely in terms of the Weyl covariant 



formalism of [55], making the Weyl transformation properties of our solution manifest. 

Having constructed these bulk metrics dual to fluid dynamics, we then proceed to 
study their causal properties. In particular, we generalize the results of || to demonstrate 
that (under appropriate assumptions) all (r = 0) singularities present in the solutions 
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presented in this paper are shielded from the boundary of AdS by a regular event horizon. 
In fact, we find a local expression for the radial location of the event horizon as a function 
of boundary fluid data, generalizing the results of J|] . We then use our construction of the 
event horizon together with the pull back of the area form on the horizon to the boundary 
(see H) to construct an entropy current for the curved space fluid dynamics constructed 
in this paper. The non-negativity of the divergence of our entropy current is guaranteed 
by the area increase theorem of general relativity. 

As an additional check and application of our results we go on to study the most general 



AdSs Kerr solutions [651 pq, p?L pq, p9L [70, 71] in fluid dynamical terms. These exact 



solutions are labelled by their mass and two angular velocities. It was already observed in 
p| that the stress tensor and thermodynamical properties of these black holes agree with 
the expectations of perfect fluid dynamics in the large mass limit. In this paper we go a 
step further; we rewrite the bulk metric of these black holes in the coordinates employed in 
our general construction of bulk duals to fluid flows. In these co-ordinates the AdSs Kerr 
solution assumes a particularly simple (and a manifestly Weyl-covariant) form 

ds 2 = -2u^dx^(dr + rA\dx X ) + r 2 'g^ lu dx tJ 'dx u 

f \ 2m ( 1 \ 1 

UyPxu^y + uj^ X uj\ u + —u^Uy \ dx^ dx u + f 1 + -^u) a i3UJ al3 J u i) u v dx ii dx v 

(1.4) 

where uj^ v is the vorticity of the fluid, T> is the fluid mechanical Weyl covariant derivative 
(associated with the Weyl connection Ax) and 1Z is the Weyl-covariantized Ricci scalar. 
In this form, this bulk metric admits a simple and explicit all orders expansion in fluid 
dynamical terms. The expansion of ( |l.4[) to second order in the derivative expansion is 
exactly reproduced by our general metric dual to fluid dynamics, upon substitution of the 
velocity field into our metric . We also verify that our formulas for the location of the 
event horizon and the local entropy current of this black hole match the exact formulas for 
the same quantities upto second order in the derivative expansion. 

As we have explained above, we have constructed gravitational solutions dual to every 
solution of the covariant Navier Stokes equations 

V M T^ = (1.5) 

on an arbitrary background manifold with metric g^ v . Let us now consider a special case 
of this equation; let 

9»u = (5oW + K) (1-6) 

for small and work to first order in h". The Navier Stokes equation may be rewritten 
as 

(v ) M (T„r = r (i.7) 



7 It is interesting to note that in the special case of solutions with vorticity lu^ v = 0, our metric exactly 
reproduces the AdS Schwarzschild blackhole. 
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where {Tq)^ u is the stress tensor of the corresponding fluid flow in the space-time with 
metric (go)uu, an d the effective forcing function f v is given by 

U = ~(T )«d a h - d.Knrh^} + ^(T )^d v h aP - VPiST&g®] 

All indices in this equation are now raised and lowered with the metric (go)^ and Vo is 
the covariant derivative with respect to the same metric. In particular, for a symmetric 
tensor Q^ v we have 



Q"J - \Q aP d v g% 



In the expression above, 5T av is the first order change in the fluid stress tensor ( 3,1S| ) 



under the variation of the metric given in ( |1.6| ), and the variation in the four- velocity 
<5u M = —^-u a u^h a j3, keeping the temperature fixed. For instance, 5T au for the perfect 
fluid stress tensor T av = T 4 (4ii Q u^ + g au ) is given by 



a pu a w 



By choosing h^ v appropriately we can produce a wide range of forcing functions 
with one qualitative restriction; is relatively mild as it is necessarily at least of first 
order in the boundary derivative expansion. Nonetheless we expect even this mild forcing 
function to be able to stir the fluid into flows with velocity differences v of unit order. The 



reason for this expectation is that that every term on the LHS of (1.7) is also of at least 
of first order in derivatives, so that factors of 1/L should 'cancel out' between the LHS 
and RHS. In below we present evidence for this assertion by presenting some simple 
steady state solutions to the equations of fluid dynamics that are pushed into flows with 
large velocities. 

The Reynolds number for flows of the conformal fluids studied in this paper is given 
by TLv; where L is the length scale of variation of velocities as above. If v can indeed be 
stirred to order unity, it follows that the Reynolds numbers of the corresponding flows are 
very large in the limit TL 3> 1 considered in this paper. Recall that fluid flows with large 
Reynolds numbers are expected to be turbulent. The discussion of this paragraph suggests 
that it should be possible to choose h^ v and correspondingly the functions j v to stir the 
boundary fluid into, for instance, steady state turbulence. The map from fluid dynamics to 
gravity, presented in this paper then yields a bulk dual to this turbulent fluid flow. It would 
be very interesting to explore this in more detail, and in particular to investigate whether 
bulk dual considerations could shed new light on the apparent universality of turbulent 
flows. 

In this introduction so far, we have described the construction of the bulk solution 
dual to the long wavelength dynamics of a field theory on an arbitrary weakly curved 
manifold. This discussion can be generalized to the bulk solution dual to an arbitrary 
fluid dynamical flow of the field theory, on an arbitrary weakly curved manifold, and with 
an arbitrary slowly varying coupling constant. In order to do this we focus on the class 
of gravitational theories which admit consistent truncation to Einstein-dilaton sector; this 
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includes all supergravity theories that result from the compactification of IIB supergravity 
down to AdSs, and so includes the duals to all well understood examples of the AdSs/CFT4 
correspondence. 

For this class of theories, we determine the dual to the the fluid dynamical motions of 
field theories with actions of the form 

S = [ y/ge~* C (1.10) 



where 0(a^) is an arbitrarily specified slowly varying function. 8 As the Lagrangian ( 1.1C| ) 



explicitly breaks translational invariance, the equations of motion obeyed by the velocity 
and temperature fields of this system are modified; in fact the Noether procedure for 



translational invariance, applied to ( 1.10 ), yields the equation 

V M T^ = e~+ C V u (f) 

T L_^_Q (L11) 

i.e. the Navier Stokes equations with an additional explicit forcing term. This equation 
also follows from a direct analysis of the Einstein-dilaton system near the boundary, upon 
employing the usual holographic formulas for T^ v and the Lagrangian (see Appendix [A]) . 

In this paper, we perturbatively determine gravity dual descriptions of fluid flows that 
obey ( 1.1 1] ) and in the process determine expressions for the stress tensor and the La- 



grangian, as a function of the velocity fields, temperature field, curvature and coupling 
constant, order by order in the derivative expansion. The procedure we employ to derive 
these results is a straightforward generalization of the procedure described above for ar- 
bitrary weakly curved boundary metric. We search for solutions of the Einstein-dilaton 
system with boundary conditions on the metric field as described above, but additionally 
require that the dilaton field <I>(x M ) asymptote to cf)(x^). The configuration ([h^) supple- 
mented with $(x M ) = 4>(x^) once again turns out to be a suitable first term in a systematic 
perturbative expansion of solutions of the Einstein-dilaton system in powers of e. This 



perturbation theory works only if the integrability condition ( 1.11 ) is obeyed. When this 
equation is met, the perturbative procedure generates expressions for T^ v and C as func- 
tions of velocities, temperatures and background fields. We have explicitly worked out this 
perturbation theory - and consequently the corresponding expressions for these quantities 
- explicitly to second order in derivatives. Under suitable assumptions regarding the late- 
time behavior of the boundary-metric and the variation of the coupling constant, we have 
also determined the location of the event horizon and an expression for the entropy current 
for this larger class of bulk configurations. 

Upon solving the gravitational equations it turns out that the expectation value of 
the Lagrangian is itself proportional to a derivative of the dilaton; as a consequence, the 



8 In sections and §^ below, we economize on space by presenting only our most general results which 
account for the effects of a varying dilaton field. The bulk metric dual to fluid flows on a curved spacetime 
with constant dilaton may easily be obtained from the results of §^ and by simply setting the boundary 
dilaton to a constant in all equations. 
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explicit forcing function in ( |1.11| ) is of second order or higher in boundary derivatives, 
and so is milder than the effective forcing function described by a metric fluctuation (|0]) . 
Nonetheless fluid flows forced by the dilaton are interesting for their bulk interpretation. 

While we have explicitly derived a duality between solution in fluid mechanics of 
strongly coupled gauge CFTs and solutions of the equations of gravity, we expect this 
duality to continue to provide a map between fluid flows in the CFT at arbitrary nonzero 
coupling and solutions to the equations of the appropriate dual classical string theory. The 
equations of fluid dynamics presented in this paper have been computed from Einstein 
gravity and so are only valid at large 'tHooft coupling, A = e^N. However, many of the 
conclusions in this paper only depend on the fact that the forcing function has a mild 
effect, if the external gravitational and dilaton fields vary slowly enough compared to the 
temperature. This feature is determined by symmetry considerations and must be valid at 
arbitrary A. 

With this in mind, consider the evolution of the hot boundary field theory with a time 
dependent dilaton that starts at a large value, dips at a particular time to a value of order 
unity, and then climbs back to a large value at later times. It seems reasonable to guess 
that the equations of fluid dynamics derived in this paper provide qualitative (though not 
quantitative) guidance to the nature of the evolution of a hot classical solution of string 
theory through regions of string scale curvature 9 . In section §0 we discuss such an evolution 
which turns out to be rather mild, suggesting the absence of qualitative surprises in at least 
this particular excursion through string scale curvatures. 

To conclude this introduction, we reemphasize that in this paper we have determined 
the gravitational dual description of the motion of a forced fluid. Forced fluid flows have 
been the subject of intense investigation for over a century |72]| and display fascinating 
dynamical behaviors all of which must be inherited by gravity. We hope that this connection 
can be used to extract useful lessons for the study of gravity and/or fluid dynamics. 

2. Long wavelength solutions of the Einstein-dilaton system in derivative 
expansion 

Consider the Einstein-dilaton system with a negative cosmological constant and boundary 
conditions as described in the introduction. In this section, we will describe how one may 
systematically solve for a class of long wavelength solutions of this system order by order 
in e = where L is the length scale of variation of the dilaton field, boundary metric, 
velocity and temperature fields and T is the local temperature of the fluid. All the results of 
this section apply to any CFT whose bulk dual description admits a consistent truncation 
to the equations of motion that follow from the Lagrangian 

S = Tok/^( fl + 12 -M (21) 

9 Recall that the curvature in string units of the corresponding string backgrounds is known to scale as 
A~4 at strong coupling and is expected to be of unit order when A is of unit order. 
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This property is true of all effective actions that are obtained from the compactification IIB 
theory on an Einstein manifold, and so of every one of the infinitely many known examples 
of the AdS/CFT correspondence. 

Notice that the Lagrangian Q2.1| ) admits a further consistent truncation to Einstein 
gravity with a negative cosmological constant if we set the dilaton field to any constant 
value. Consequently, when our boundary conditions allow us to set the dilaton to a constant 
the equations from (2.1) further simplify. All the results of this section and the next apply 



with greater universality in this limit - they apply to any CFT whose dual description is a 
two derivative theory of gravity interacting with other fields. 

2.1 Uniform Branes in the Einstein-dilaton System 

The equations of motion that follow from Lagrangian (2.1) are given by 



_ R . Of - 20, R MN - + iGm _ „ < 2 - 2 > 
V 2 $ = 0. 

There exists a well known four parameter class of exact solutions to these equations of 
motion 

ds 2 = -2u^dx^dr - r 2 f{br)u^u v dx^dx u + r 2 V '^dx^ dx u 
$ = $o 

T> J- ( 2 - 3 ) 

where is an arbitrary constant matrix of signature (—1, 1, 1,1), b is a constant and 
is a constant unit normalized velocity vector : u^u u g^ u = — 1. 

2.1.1 Regulation and Weyl Frames 

In this subsection, we pause to recall a well known subtlety in the boundary interpretation 
(via the AdS/CFT correspondence) of locally asymptotically AdS metrics like (|2.3|). This 
interpretative subtlety will play no role in the calculation we describe in this section or the 
next. However, it will permit a stringent test on the self-consistency of our final results in 

We begin by noting that the metric ( [2.3| ) asymptotes to 

ds 2 = -2u^dx^dr + r 2 g^dx^dx" (2.4) 

and according to the usual rules of the AdS/CFT correspondence, describes a state in the 
CFT formulated on a space with any of the infinite numbers of metrics that are Weyl 
equivalent to g^ v . . Next, we note that in order to find the dual interpretation of any 
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asymptotically AdS space, it is convenient to regulate the solution near its boundary and 
the bulk solution may then be regarded as a state on the field theory on a base manifold 
whose metric is proportional to the induced metric on this regulated boundary. 

In particular, we may choose to regulate the boundary of (2.4) on slices of large 



constant r. With this choice the dual CFT resides on a space whose metric may be taken 
to be precisely g^, and ( |2.3[ ) in fact represents the CFT on a space with metric g^ u uniform 
temperature T = —r and in uniform motion at velocity u^. Alternatively, we may choose 
to regulate the boundary of ( |2.3D along slices of constant f = e~ x ^ x ^r for any arbitrary 
function x{ x>1 )- I n order to interpret this new slicing it is convenient to rewrite the metric 
in terms of f; Asymptotically, we have 

ds 2 = -lu^df + r 2 ~g liv dx il dx v (2.5) 

where u M = e x u^, and g^ v = e 2x g^ u . It follows that the constant r slicing of Q2.3| ) describes 
a state the CFT on with spacetime varying background metric g^, and with spacetime 
dependent velocities and temperatures T = i = eT x T for arbitrary xi xf1 )- 

The fact that different slicings of the same supergravity solution may in fact be inter- 
preted as states of the same theory in distinct though Weyl equivalent background metrics, 
reflects the conformal invariance of the dual field theory. As we have seen above, the 
temperature and velocity fields T, and transform homogeneously under Weyl transfor- 
mations, with weights +1 and —1 respectively. 

2.2 Long Wavelength solutions with slowly varying background fields 

In this subsection, we will explain the procedure that we use to construct a class of solutions 
to the Einstein-dilaton action ( |2.lD that asymptote at large r to the metric and dilaton 

ds 2 = -lu^x^dx^dr + r 2 g iiv {x^)dx^dx 1 ' 

$ = 4>{x) 

for arbitrary long wavelength functions 4>{x^) and g fiu (x^). The solutions we construct 
are dual to fluid dynamics in a field theory with an arbitrary long wavelength background 
metric and coupling constant given by g^ix^) and 0(x^). The method we employ in this 
paper closely follows the method in 

2.2.1 The Zeroth Order ansatz 

We first note that the ansatz 

ds 2 = -2u^)dx»dr + r 2 f{b{x tJ ')r)u^u v dx ti dx v + P^dx^dx" 
$ = ftx") 

?V = gA^) + *v(^KO M ) (2 ' 7) 

for arbitrary functions u^{x^) and b(x^) has several attractive features. First, it asymptotes 
at large r to (|2.5| ). Second, all singularities in ( |2.7D (which occur at r = 0) are shielded 
from the boundary by a regular event horizon. Third, this configuration solves the Einstein- 
dilaton equations of motion for constant g^ u , (j), and b. 
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2.2.2 Setting up a derivative expansion 

As in Q , we will use fl2,3| ) as the first term in a systematic derivative expansion of a solution 
to Einstein's equations. The small parameter that justifies this procedure is the inverse 
length scale of variation of the functions g^, <j>, and b (each of which is assumed to vary 
on the same scale) times the function b. Denoting this parameter by e (see fl|] for more 
details) we plug the expansion 



V u ~ / j \w 

m=0 /„ 
oo \ Z - 



m=0 



into the Einstein-dilaton system and solve Einstein's equations perturbatively in e. As we 
will see, the expansion of these equations to m th order in e will allow us to determine the 
functions G™ and $ m . 

2.2.3 Gauge Choice and Physical Interpretation 

Following Q , in the bulk of this paper we work with the gauge choice 

G rr = 0, G ril oc U/M Tr ((G°) _1 G m ) =0 (m > 0) (2.9) 

The first two of these gauge conditions have been physically interpreted in ||. As ex- 
plained in that paper, lines of constant x M are geodesies in any metric subject to this gauge 
condition. The last gauge condition in ( |2.9D was chosen arbitrarily and turns out to have 
no interesting geometrical consequence. However, it is possible to choose a more natural 
gauge condition with a good geometric interpretation - 

G rr = 0, G rfl = (2.10) 

(2.10) continues to ensure that lines of constant are geodesies but also ensures that the 



coordinate r is an affine parameter along these geodesies. 

In order to permit easy comparison with the results of Q, we will work with gauge 
( |2.9| )the bulk of this paper. However, we will also briefly indicate how our final results may 
be transformed into ( |2.10| ), which will prove useful when comparing with explicit black 
hole solutions. 



2.2.4 Constraint Equations 

As in W it turns out that only 15 of the 16 Einstein plus dilaton equations are actually 
independent. These 15 equations may be separated into 4 constraint equations (see |TJ) 
and 11 dynamical equations. 

In Appendix |A|, we show from a direct analysis of the Einstein-dilaton system near the 
boundary that the four constraint equations take the form 



(2.11) 



where the stress tensor T^ u and Lagrangian C are defined by in terms of gravity data by 
the usual formulas of AdS/CFT 10 

ltorG 5 T£ = lim r 4 (2(^/1^ - K£) 

-fty-\ (v"*V„* - f (V$) 2 )) (2.13) 

16vrG 5 e-* C = - lim r 4 ( <9 n $ + ^V 2 $ 

where n M is its outward pointing unit normal of the regulated boundary and h^ u is its 
induced metric which leads to the covariant derivative V and the corresponding Einstein 
tensor Qy. The extrinsic curvature of the regulated boundary is defined via the normal 
lie-derivative of the induced metric - = \Z n h ilv and d n is the partial derivative along 
. All the indices in the above formulas are raised using the induced metric. 
As each term in these equations has an explicit boundary derivative, the constraint 
equations are special from the viewpoint of the boundary derivative expansion. This is 
because each boundary derivative pulls down an additional power of e; consequently in the 
expansion of the constraint equations to order m we find contributions from the functions 
G™ v and <& m only for n < (m — 1). Consequently, the constraint equations at m th order 
do not aid in determining the unknown functions G™ u and <£ m ; they instead impose a 
constraint on the solution obtained upto (m — l) th order. As this solution has already 
been determined (by the perturbation theory to (m — l) th order) in terms of u^{x^) and 
T(x^), these equations effectively reduce to equations of motion for these fluid dynamical 
fields. In other words, at any given order in perturbation theory, the constraint equations 
are simply the equations of fluid dynamics. 

2.2.5 Dynamical Equations 

The expansion of the remaining 11 independent dynamical equations, to order e m , yields 
equations that may be used to determine the unknown functions G™, <& m . 

The nature of the resultant equations is described in pj. These equations are ultralocal 
in the variables (they only contain derivatives of r). Consequently they may be solved 
independently point by point in and are linear ordinary differential equations in the 
variable r at each x^. At any given x^ these equations take the form 

N($ m ) = s m 

where the homogeneous operators M and iV are independent of m (they are the same at 
every order in perturbation theory). 



10 Note that the counterterm subtractions to the stress tensor and Lagrangian respectively in ( 2.13 ) are 
simply the derivative with respect to g^v and <j> of the counterterm action |73| 



S ctr = / Vh - - -{d4>Y (2.12) 
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2.2.6 The Homogeneous Operators 

Let us describe the differential operators M and N in more detail. Focus on the equations 



( 2.14 ) at a particular field theory point y^. It turns out that the operators M and N 
are simply the equation for linearized radial fluctuations of the metric and dilaton in the 
background of the uniform brane ( [2.3; ) with constant boundary metric guv(y^), constant 
dilaton (f>(y^), constant velocity u^(y^) and constant inverse temperature b(y^). In other 
words the homogeneous operators M and N are simply the operators that act on linearized 
radial fluctuations of the metric about the background ( |2.3[) . 

We can now perform a local co-ordinate transformation to go to Riemann normal co- 
ordinates where the metric g^ v = rj^ . We still have some more co-ordinate freedom param- 
eterized by the set of Lorentz transformations which preserve this form of the metric. These 
boost transformations may be used to set = (1,0,0,0). Finally, a scale transformation 
— > Ax^ and r — > j- may be used to set b to unity. Making this choice of coordinates 
for any given point y^ we see background metric ( |2.3| ) turns into the metric of a uniform 
black brane at rest with constant temperature T = i with the usual flat boundary metric. 
However this was precisely the background metric encountered in the perturbative proce- 
dure described in [Q]. It follows immediately that operator M is identical to the operator 
H in equation 3.4 of Q . 

In the same coordinates the operator N is easily determined; it is simply given by 



N = \d r (r\l - \)d r ) , (2.15) 



5 ( 

the radial part of the minimally coupled equation in a black brane at b = 1. 

As was explained in the equation MG^ = s^ v may be solved by integration for 
an arbitrary source s^ u . As we will explain below, the same is true for the operator N. 
Consequently dynamical equations may be solved by integration at each order. As in (l|, 
in this paper we will solve these equations subject to two boundary conditions 

1. G™ v and <£ m are well behaved (analytic) away from the r = singularity. 

2. $ m and r 2 G^, each go to zero as r — ► oo. 

The first condition requires no explanation, while the second one ensures that our correc- 
tions to the metric and dilaton fields do not alter the asymptotic form ( |2.6[ ) of our full 
series solution. 



2.2.7 Source Terms 

We now turn to a brief discussion of the source terms s™ u and s m . As in Q each of these 
source terms is determined in terms of lower order solutions in perturbation theory, and so 
may be expressed as local functions of g^ tv (x),(j)(x),u tM (x) and b(x). s™ and s m are each 
of m th order in derivatives of these quantities. The dependence of s™ on derivatives of 
b and was already described (and determined to second order) in ||]. Below we will 
also determine the dependence of s™ on derivatives of the metric and dilaton (for m < 2) 
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and compute s m for m < 2. In this subsubsection we explain the general structure of our 
results. 

Let us first study the dependence of source terms at on derivatives of the metric, 
at the same point, for m < 2. For the purpose of this calculation we find it useful to 
work with Riemann normal coordinates centered about the point y^ in the boundary field 
theory directions. In these coordinates the expansion of the boundary metric g^ u about y^ 
starts at second order. As a consequence s\ v and s 1 each are independent of derivatives of 
9[iu- While s 2 ^ receives contributions proportional to the boundary curvature tensor, s 2 
(like all dilaton source terms) is necessarily proportional to a derivative of some order of 
the dilaton field 11 Consequently s 3 is the first dilaton source term that receives curvature 
dependent contributions; s 2 is independent of boundary curvatures. 

We now turn to the dependence of s™ v and s m on derivatives of the boundary value 
of the dilaton for m < 2. Let us first focus on s m . As we have argued each of s m must 
be proportional to at least one derivative of the boundary dilaton. From symmetry it 
then follows that s 1 oc u^dcf). On the other hand s 2 can and does receive contributions 
proportional to both vP"u h 'd^dycj) as well as terms proportional to one derivative of a velocity 
contracted with a derivative of 4>. 

We now discuss the contribution of derivatives of the boundary value of the dilaton, 
4>(x) to s™ u . It follows from (|2.2|), the dilaton source term for the metric equation is propor- 
tional to two derivatives of the dilaton. As the dilaton is constant (has no r dependence) in 
the uniform brane solution, and (as we have seen in the previous paragraph) is proportional 
to a derivative of the boundary value of <j) at first order, it follows that dilaton contributions 
to vanishes, and the contributions of the dilaton to s 2 ^ are schematically proportional 
to d^cj)d v (j). 

2.2.8 The Stress Tensor, Lagrangian and Fluid Equations of motion. 

As we have explained above, the dynamical equations upto m th order in e may be used to 
solve for the metric and dilaton, and consequently the stress tensor and the Lagrangian 
(from ( |2,13| )) to the same order. By plugging into ( 2,11| ), this information also deter- 



mines the fluid dynamical equation (including the forcing functions) with (m + 1) or fewer 
derivatives. 

This concludes our brief review of the general structure of the long wavelength pertur- 
bation theory we perform in this paper. We refer the reader to [Q] for a fuller description 
of the procedure. 



3. Explicit Results Upto Second Order 

In this section we present explicit formulas for the metric and dilaton, as a function of 
boundary metric, dilaton, velocity and temperature fields, to second order in the boundary 
derivative expansion. We have obtained these results by implementing the perturbative 
procedure described in the previous section and in [Q]. 

11 This follows from the fact that Einstein gravity is a consistent truncation of the Einstein-dilaton system. 
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3.1 The metric and dilaton at first order 

As we have described in the previous section, the metric dual to fluid dynamics at first order 
is already completely determined by the results of Q. s^^y 11 ) receives no contribution from 
first derivatives of the boundary metric or boundary dilaton at y^, and the metric G 1 ^ is 
simply by the naive boundary covariantization of equation of equation 4.24 of Q. 

The dilaton field $ x is nonzero and requires a new - though very simple - calculation 
to determine. 

Using coordinates in which g^uiv^) = Vnui = (—1,1,1,1) and b(y^) = 1, the 
equation for <I> 1 (y M ) is 

8 r (V 5 (l - ^)d r ^A + d r (r 3 d v 0) = (3.1) 

This equation may be integrated trivially. The arbitrary solution to this differential equa- 
tion is given by 

f°° 1 r°° r 3 — 1 

* W, m) +9A L (3 - 2) 

Our boundary condition at infinity forces c\ = while the requirement of regularity at the 
horizon r = 1 sets C2 = 12 . Consequently we conclude 

f°° r 3 — 1 f°° r 3 — 1 

* 1=8 W r ^wr-W, ?w) (3 - 3) 

3.2 Solution at second order 

The computations required to determine the metric and dilaton field at second order, while 
involved in practice, are a straightforward generalization of the calculations presented in (l| . 
We have performed these calculations with the aid of the symbolic manipulation program 
Mathematica. In this section we simply record our final results. 
The metric upto second order is given by 

ds 2 = - 2u fl (x f *)dx fl dr + r 2 f{b{x^)r)u^u u dx^dx v + V^dx^dx" 
2 b r 2 F(br)a^ u + 9 u^u u - r(a^u u + a u u^j dx^dx v 
+3 b 2 H u^dx^dr ^ ^ 

+ ^r 2 b 2 H V^u + -^2 K u n u v + ^2 (JyVv + Jf u n) + r 2 b 2 a^ dx^dx u 

$ = + u.d<p r dx^-\ + Sfhf{br) + Sfhf{br) 

J rb x J\ x ) 

In this equation the first line is simply the ansatz (2.7). The second line records corrections 
to this metric at first order in e, while the third and the fourth lines record the second 
order corrections to this ansatz. 



12 Note that, for any non-zero value of C2, the dilaton has a logarithmic divergence at r = 1. This (co- 
ordinate invariant) singularity is not shielded by any horizon and hence we will consider it to be physically 
unacceptable. We will show in the later sections that when ci = 0, our solutions have regular event horizon 
that passes through the neighbourhood of r = 1. 
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(3.4) - the metric dual to fluid flows upto second order in the derivative expansion - 
is one of the main results of this paper. In the rest of this section we will systematically 
define all the previously undefined functions functions that appear in ( |3.4|) . We will start 
by defining all scalar functions of the radial coordinate r that appear in (2.7), and then 
turn to the definition of the index valued forms that these functions multiply. 

A3 



The only undefined function of r in the second line of (|3,4| ) is F(r) which is given by 



F(r) 



In 



l+r) 2 (l+r 2 ) 



2 arctan(r) + tt 



(3.5) 



The undefined functions on the third fourth and fifth line of the same equation are defined 

as 



H = /i (1) (6r) 64 + /i (2) (6r) 65 + /i (3) (6r) Si* 

K = k^\br) 64 + k^(br) 65 + k^\br) S + k^(br)&i* + k (5 \br)Gz 
J M = j«(&r) B- + f 2 \br) B^+j^{br)Bi 



a> = ai{br)T^ u + a 5 (6r) (T 5 ) Mi , 

+ a 6 (br) (Tq)^ + a 7 (br) (T 7 )^ v + a 8 (6r) + a 9 (6r)(T^) ^ 



(3.6) 



where 



h<®(r 



1 



12r 2 
1 



6r 2 



+ 



00 dx 
x 5 



dy y 4 -W h ( y ) 



r 4 + 3 r 4 + 3 



+ 



96r 4 
r 4 + 3 
48r 4 



TT 



tan l (r) lnr 



48r 4 



ln(l + r) + 



1 



32r 4 



12 



ln(l+r 2 ) + 



2 



-2r 3 + 3r 2 + 1 
48T 4 



(3.7) 



13 As explained in jjj], we obtain -F(r) by integrating a second-order differential equation. One of the 
integration constants of this equation are fixed the requirement of regularity of F(r) at r — 1. For a general 
value of this constant, F(r) has a logarithmic singularity at r = 1. This singularity is physical rather than 
a coordinate artifact, and is hence unacceptable. We have checked the last statement by computing the 
curvature invariant Rmn R mn ; this quantity in general has a pole type singularity, proportional to a^ lv a l ' tv 
at r = 1. 
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r 2 r°° ( 

k<U(r) = -— - J \ 12x 3 h^\x) + (3x 4 - 1) 
fc ( 2 )( r ) = _ J°° (i2xWXx) + (3x 4 - l) dk{2){x 



dhP-\x) l + 2x 4 x 
dx 6x 3 6 



fc(3)( r ) =r 2/ 2 

+ 



32r 4 
r 8_j 



1 r 

7T 



1 tan _1 (r) + 7 (l - r 4 ) lnr 



16r 4 

Ml + r) + ^ Ml + r2 ) 



+ 

jfe(5)( r ) =r 2/12 



16r 4 ' '' ' 32r 4 

6r 7 + 5r 6 + 3r 4 + 2r 3 _ 3r 2 _ 1 



48r 4 



r 2 />oo /-oo / 

J<"(r)- s -jf ■bx'jf *( 

./r 



18y 3 (y + l)(y 2 + l) 9y 3 



,-(2) (r) 

ai(r) = 
05 (r) = 
a 6 (r) = 

fl7(r) = l/ 



poo />oo 
„3 



dxx^ ^Vl«'/*0/ + J)(// 2 + D 
1-r 4 



f 

Jr 

f 

Jr 



16 

dx 



it - 2tarT 1 (r) + 21n(l + r) - ln(l + r 2 )) 

My) + 11 



r 2 (l + 2r) 



a 8 (r) = - jf 
a 9 (r 



x(x 4 - 


1) 


dx 




x(x 4 - 


1) 


dx 




x(x 4 — 


1) 


(ix 




x(x 4 — 


1) 


dx 





dy 2y 



3yF(y) 



p{y) + 5 

^ dy2y^2 



p(y) + 5 



GyF(y) 



p(y) + 1 

p(y) + 5_ 



6yF(y) 



x(x 4 

Mi + 



2y 



W ft (r) = 
W fc (r) = 



4 (r 2 + r + l) 2 - 2(3r 2 + 2r + l)F(r) 
3 r(r + l) 2 (r 2 + l) 2 

2 4(r 2 + r + l)(3r 4 - l)F(r) - (2r 5 + 2r 4 + 2r 3 - r - 1) 
3 



r(r + l)(r 2 + 1) 



h 



p(r) = 2r 3 + 2r 2 + 2r - 3 
dx 



(i) 



+ |.tan- 1 (y) + 



y(l + y)(l + y 2 )-2y 



f°° dx f 
(r) = X x(x 4 - 1) A 



3(l + y)(l + y 2 ) 
3 



(3.8) 



(3.9) 



(3.10) 



(3.11) 
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We now turn to denning all the terms that carry boundary index structure in (|2,7|). 
These terms are all expressed in terms of fixed numbers of boundary derivatives of the 
velocity, metric and boundary dilaton fields. 

Terms with a single boundary derivative 



6 = V a u a 

<V = \v» a V uP (V a up + Vpu a ) - ^P^O 

The quantity defined here does not appear in the first order correction to the ansatz 
metric, but does appear, multiplied by other first order terms appears in the second order 
metric. 

We now describe all terms with two boundary derivatives. Following we sub-classify 
these terms as scalar like, vector like or tensor like, depending on their transformation 
properties under the 50(3) rotation group that is left unbroken by the velocity (see Q 
for more details) 

Scalar terms with two derivatives 

-t(s 3 -Ri) + 2 ©i-^©3 

62 = 



64 = V" 

©5 = <v °^ 



(3.13) 



2 



Mi 



where 



Si = a^ 

& 3 = e 2 

s 3 = \v«N a Vp b ( 3 - 14 ) 



Vector terms with two derivatives 



B°° = 4 (10 V4 + v 5 + 32Ji - 35J2 - 6533 + 9 R 2 ) 

B fin = 9 (20 v 4 - 5 2J2 - 6 2J3 + 20 R 2 ) (3.15) 



b* = u.vcp pjy t 
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where 



(V5)v 



9 
5 

a a au 
1 

~2 



(3.16) 



^ tit* Rjxoryf} 



Tensor terms with two derivatives 



(^5)^ 
(^6) iiv 



+ 3(2:4)^ + (1 3 ) 



1 



TV ®5 



K V » ( R "P + 2u7 ^7^/3) - (Rap + n A i? 7QA/3 ) 



(3.17) 



where 



a,, a 



o?V ©I 



(I3W = ( u - v ) ( v « n /3 + v /^«) ~ l^r af3 (u.v)(v a u p ) ( 3 - 18 ) 



Here V denotes the covariant derivative in the curved boundary metric, R a /3-yu is the 
Riemann tensor and R^ is the Ricci tensor. 

The metric and dilaton configuration presented in this subsection solves the coupled 
Einstein-dilaton equations of motion provided that the temperature and velocity fields obey 
the equation of motion presented in the next subsection. 



3.3 Shift of Gauge 

While we have (for historical reasons) presented our final metric gauge choice ( |2.9[ ), it 
would have been more natural, in some respects, to work with the gauge choice (|2.10| ). 
The variable change that converts from ( |2.9[) to ( |2.10| ) is simply given, to second order, by 
df = dr(l — ^b 2 H(br)). This variable change may be integrated. The metric expressed in 
terms of the variable f will be in the gauge ( |2.10| ) at second order. 



- 19 - 



3.4 The Stress Tensor, Lagrangian and Forcing function at second order 

The stress tensor and expectation value of the 'Lagrangian' may be determined for the field 
theory configurations dual to the solutions of the previous subsection using the formulas 
(|2l3|) . We find 



lfaGsT^ = (tt T) 4 + 4 u^u v ) - 2 (vr T) 3 a*" 

+ (ttT) 2 ((^) (T 7 )^ + 2(T 6 )^ + (2-ln2)T^ + C^-ir^ 

-l*KG s e-*£ = (7tT) 3 u.d<P + (vrT) 2 Q«sf + 

(3.19) 

Further, the spacetime configuration presented in the previous section is a solution 
to the equations of motion if and only if the velocity and temperature fields obey the 
constraint 

V, t T^ = f v 

a. (3-20) 



4. Constraints from Weyl Covariance 



In the previous section we have determined the bulk metric dual to a particular fluid flow 
Ufj,(x) on a boundary manifold with a particular background metric g l _ iv {x^) and boundary 
dilaton field (f>(x^). However, as we have explained in section [2|, this must be the same 
as the bulk geometry dual to the velocity field u^(x^) = e x u fl (x 11 ) on the boundary space 
with metric g^ix^) = e 2x g^ v (x^) and the boundary dilaton unchanged. This equality of 
spaces, which is not algebraically automatic in the formulas of our perturbation theory, 
constitutes a tight algebraic check on our procedure. 

In this section we will demonstrate that the metric described in the previous section 
passes this test. In more detail we will demonstrate that the metrics dual to fluid dynamics 
must are invariant under the simultaneous replacements 14 



u^^e x u^, g^tu^e x g flu , b -> e x b, r -> e x f, cf> -y <j> (4.1) 

This property of our solutions in particular implies the Weyl covariance of the fluid 
dynamical stress tensor (and Lagrangian) that follows from our solutions. 



14 It may be useful to reiterate the logic that underlies this test. Let us imagine we have solved the problem 
described in the previous subsection for a given background metric and velocity and temperature fields 
Ufj.{x) and b{x). Upon performing the coordinate transformation r = re - *^' (as described above) for a 



slowly varying function x( x ) we have a gravitational background of the form (2.7) with the new metric 
velocity and temperature functions w M = e x u M and — e 2x g MI , and b = e x b. However we could also have 
directly solved for the metric dual to fluid dynamics with this data. 



- 20 - 



4.1 Weyl Covariant Formalism and Independent Weyl Covariant Tensors 



Note that the transformation (|4.l| ) is simply a Weyl transformation in the boundary field 
theory directions. Consequently, in order to investigate the invariance properties of our 
solutions under the transformation (f4.1[ ) we find it convenient to employ the manifestly 
Weyl-covariant formalism for hydrodynamics that was developed in f55fl. The main tech- 



nical innovation of [55] was the introduction of a Weyl-covariant derivative, whose action 



on an arbitrary tensor field is defined by 

+ [gxaA" - 5%A a - 6£A X ] QZ'.W + ■■■ (4.2) 
- [gxuA a - 5^A U - % Ax] Q&: - • • • 

where the Weyl-connection An is related to the fluid velocity via the relation 



A /1 = a fl - -j—^Uti ( 4 - 3 ) 



As a technical prerequisite to our main goal, in the rest of this section we employ 
this Weyl covariant formalism to list the most general Weyl invariant scalars, transverse 
vectors and symmetric traceless transverse tensors in hydrodynamics that involve no more 
than second order derivatives. 15 We perform this listing, taking particular care to account 
for the equations of motion. In other words tensor fields that are equivalent on-shell are 
counted only once in our listing. 

Let us first start by eliminating an easily dealt with redundancy. As an arbitrary 
function of the dilaton is Weyl- invariant, it follows that we can generate new Weyl invariants 
by multiplying old ones by functions of <f> to get new Weyl-invariant observables. We shall 
omit the observables formed this way from the lists below, but take care to account for 
them later. 

To begin, let us start with the basic quantities of hydrodynamics - the fluid temperature 
T and the fluid velocity u^ 16 . The former is a Weyl-covariant scalar with conformal weight 
unity and the latter is a Weyl-covariant vector with conformal weight unity. It follows that, 
at the zero derivative level, there are no non-trivial Weyl-invariant scalars, no transverse 
vector or symmetric traceless transverse tensors. 

17 



15 We will restrict attention to fluid dynamics in 3 + 1 dimensions. 

16 We will assume that there are no other conserved charges except the energy momentum tensor. 
17 We shall follow the notations of ^| in the rest of this section (except for the curvature tensors which 
differ by a sign from the curvature tensors in [55[ |. In particular, we recall the following definitions 

K = R + 6X7 X A X - 6^ A ^ A ; 2V*„ = <x M „ + ui^ ( 4 - 4 ) 



Note that in a flat spacetime, R is zero but 1Z is not. 
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At one derivative level T^u^T)^ is the only Weyl invariant scalar. The only Weyl 
invariant transverse vector is P^T> u <p. Finally, the only Weyl-invariant transverse pseudo- 
vector and only one Weyl-invariant symmetric traceless transverse tensor Ta^ u . 

At the two derivative level, there are seven independent Weyl-invariant scalars 

T- 2 P^ v V^V u <p, T- 2 u^u v V^V u (j), T- 2 P^V^V u (j) and T- 2 u^u v V^V v 4> 

(4.5) 

one Weyl-invariant pseudo-scalar T^ 2 l^T>^cj) and six independent Weyl-invariant transverse 
vectors 

T- l P;V x a u x , T- l P»V x u u x , T- l Plu x V v V x 4>, T^P^V^V^, 



T- l a^V x <P, and T^cj^V^ 



and two Weyl-invariant transverse pseudo-vectors T l u x T>\l il and T 1 l ll u x T>x4> . 
There are eight Weyl-invariant symmetric traceless transverse tensors - 

u A £> A o>, (j^ u u x Vx4), C mv pu a u^ ', uj^axu + ^u X cr\ fl , 

1 " 



2 



o 

pa pP , pa p/3 pa/3 p 



2 

pa p/3 i pa p/3 pa/3 p 



\ Pfiu a , A , Puv a/3 

On <?\v y o a pa p and w M w Ai , + -y- w Q /3^ p . 



(4.6) 



and five Weyl-invariant symmetric traceless transverse pseudo-tensors 18 

2P 

e al3X »C a p ua u°ux + e a ? x v C a ^ a u°ux, e^V^W a Xu + ^\u a V p <f>a^ 

and i e a/3 % e^C^ u A « CT . 

4.2 Manifestly Weyl covariant form of the fluid dynamical metric 

In either of the gauges ( |2.9| ) or ( |2.10| ) employed in this paper, a bulk derivative metric that 
is invariant under the transformations ([0]) must be constructed out of 

1. The boundary Weyl invariant tensors listed in the previous subsection. 

2. The Weyl-invariant vector-valued 1— form u fl (dr + rAxdx x ) 

3. Arbitrary functions of the Weyl-invariant scalar £ = r/{irT) 



18 Note that the last tensor is not really a pseudotensor - but we count it along with the other tensors 
which need e a/3Xl1 for their definition. 
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. Consequently the metric must take the form 

ds 2 = -2Wi u^dx^{dr + rA x dx x ) 

+ [r 2 {W 2 g^u + W&ipUy) + r {W 4fl u u + u^W Au ) + W 5ia/ ] dx^dx". 

where W%, W% and W% are Weyl-invariant scalars, W^n is a Weyl-invariant transverse vector 
and W§uv is a Weyl-invariant transverse traceless tensor. It follows that the W%s can be 
expressed as functions of £ and the Weyl covariant observables in fluid dynamics with 
appropriate weights. 

In the rest of this section we will check that our metric (|3.4| ) is indeed of this form. 
We will also compute all the functions that enter into (|4.7| ) 

4.3 Weyl Covariant form of the second order metric 

In order to cast the metric in this paper in the form given by ( |4.7| ) , we have to first rewrite 
the quantities appearing in the metric in a Weyl-covariant form. We have the following 
relations in the flat spacetime which identify the Weyl-covariant forms appearing in the 
second-order metric - 



; Q/3 ; 65 = cj a ^- sf = iu»u v v»v v fr sf = i 

4, 2^ 2 - 2 1. 



3 -(s 3 - Ri) + 2©i - ^63 = ^a af sa a(3 - -oo a ^ - ^Jl 



3 p 12 p 6 



(4- 



These can be used to obtain 



B« n = 5AP»V x a x + 90P»V x u x (4.9) 
B^ = P»u x V u 0V x 

Hence, all the second-order scalar and the vector contributions to the metric can be written 
in terms of three Weyl-covariant scalars a a pa a/3 ,oj a pu! a/3 and 1Z and two Weyl-covariant 
vectors £>a°> A an d V x lo^ x . 

Further, we will enumerate below the covariant forms which appear in the tensor 
contributions to the metric 

X M „ = u x V x a^ v 

pfiu 

= 4(uAu/" + — <A Q/3 ) 

= ^ x a Xv - —o^o al3 (4.10) 
T 7M! ,= -2{^ x a Xu +uo\a X ^) 
(T^ = P«pPv a <t> V p cf> - \p aP P^V a (t> Vrf 



- 23 - 



Comparing it with the metric form in 4.7 , we find that the metric can be written in 
the form 



ds 2 = -2Wi u^dx^(dr + rA\dx x ) 

+ [r 2 (W 2 g„u + W 3 u M u„) + r (W^u u + Ufl W 4u ) + W^ v ] dx»dx v . 



where 



Wi = l + 



+ . 



Wo 



1 F 1 {Qcj^ + 



W 3 = £ 



6r 2 
6r 2 



+ ... 



w 



4ju 



£2 r 

F 3 (Qp»v x a u x + p;iW + 2r 2 J (3) (0^> A ^^A^ 

2r 



+ • 



= 2r£F(£)<V + W 52 i(£K2?a<V + ^WO K V ^ + ^ ^ 



+ W 523 (0 ( o x o Xv 



P 



+ -y- w Q/3 u; p 



+ e 2 in(i + r 2 ) 



C^pvTvP + ^ I P«PP - -P ap P^ ) V a <p Vrf 



1 



jet p/3 



1 



up-. 



+ 



(4.11) 



(4.12) 



with 



no = 7 



In 



2 arctan(£) + tt 



, , p2 f°° dx [°° ydy [1 + 2y(l + y + y 2 ) + 2y 2 (l + 2y + 3y 2 )F(?/)] 

lUj 5 A (y + l) 2 (y 2 + l) 2 

1 f°° 4(x 2 + a; + l)(3x 4 - l)F(x) - 2x 3 {x 2 + x + l) + x + l 



F 2 (0 = 1 



20 



x(x + l)(x 2 + 1) 



4 r°° r c 
Fs(0 = 1 + ^ J x dx j 



dy 



y3( y + l)( y 2 + l) 



(4.13) 



and 



^521 
W522 
W523 



dx 


x(x 4 — 


1) 


dx 




x(x 4 — 


1) 


dx 




x(x 4 — 


1) 



dy 2y 
dy 2y 
dy 4y 



1 + 3y + 3y 2 + 3y j 

(v + W + D - 3yF(y) 
-l + y + y 2 + y 3 

-1-?/ + 2?/ 3 + 2?/ 4 + 2t/ 5 
y 2 (y + l)(y 2 + l) 



3yF(y) 



(4.14) 



The functions and 

are as given m \p. l\). 
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Further, the dilaton field can be written in terms of the boundary values as 
$ = <P(x) + u^V^j™ dx^j^ + ShfiftuWVfPvt + hfi^V^V^ (4.15) 

where the functions K\ (£) are as given in (|3.7|). 



4.4 Constraints on the Stress Tensor and Lagrangian from Weyl Covariance 

In ( |3.19| ) we have presented formulas for the stress tensor and the expectation value of the 
dilaton upto second order in the boundary derivative expansion. In the notation of this 
section 

16ttG 5 T^ = (tt T) 4 (g^ + 4 u»u v ) - 2 (tt T) 3 g» v + (2 - In 2) (ttT) 2 u x V x g» v 
- ln2(^T) 2 {^\a Xv + u> u \<r x ^) 

+ 2(vrT) VV A " - —a^a^ + C^u a up) 

-l6iTG 5 e^C = (ttT) 3 u.V(j) + (V 2 cp + In 2 (u.V) 2 <j>) 

Aspects of this result could have been predicted immediately from the requirement of 
Weyl covariance. In particular it follows immediately from the listing of Weyl invariants 
in subsection 4J. that at the one derivative level 19 the stress tensor and the Lagrangian 
are proportional respectively to (tt T) 3 and (ttT) 3 u.d<f> respectively. Of course the 



particular coefficients in ( 4.16 ) required a calculation to determine. In a similar manner 
the two derivative contributions to the stress tensor and Lagrangian are constrained by 
Weyl invariance to be linear combinations of the forms listed in ( |4.6| ) and ( [4.5[ ) respec- 
tively. Again the particular coefficients in ( |4,16| ) require knowledge of the bulk equation 
of motion, and follow from explicit computation. Note in particular that the Lagrangian 
in ( 4.16| ), when expressed as a linear combination of the terms in Q4.5| ), turns out to 
have no terms proportional to the Weyl invariant scalars T~ 2 a^ u a^ u , T~ 2 uj^ u uj^ u ', T~ 2 1Z , 
T~ 2 P tJ-u T> fl (j)'D v (J) and T~ 2 u IJ, u v T> il <j)'D u (j). Similarly, the stress tensor, when rewritten as a 



linear combination of the terms in (|4.6|), has no terms proportional to any of an V u x T> x 



pa pP i pa pP 2 pa/3 p 

r a r v ~r r v it fx 3 _r r ^ v 



T> a T>p(j), or uj^ujxv -\ — (p u a pu a P. Above we have already 
remarked that general structural features of the Einstein-dilaton system (e.g. the possi- 
bility of a consistent truncation to pure gravity) explain several of these terms from the 
expressions for the bulk metric and hence Lagrangian and the stress tensor. 



5. Causal Structure and Entropy Current 

In this brief section we will generalize the construction of || to determine the location 
of the event horizon of the metric ( |3.4D (under certain assumptions about the late time 

19 Interestingly, the fact that the Lagrangian vanishes at the zero derivative level is not automatic from 
Weyl invariance, but follows instead from the additional dynamical input that Einstein gravity is a consistent 
truncation of the Einstein-dilaton system. Once we take a' corrections into account, we do expect a 
contribution to the Lagrangian that is simply proportional to T 4 . 
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behavior of the functions in this metric) and also to derive an explicit expression for a 
positive divergence entropy current for the dual fluid flow. 

For the purposes of this section we restrict our attention to solutions of fluid dynamics 
for which the late time metric, dilaton, velocity and temperature fields all settle down to 
time independent constant values; we further assume that the the event horizon of this 
asymptotic late time spacetime is given by r = l/b(x^). The last assumption certainly 
true for solutions that approach black branes in flat space or rotating black holes in global 
AdS space at late times. Under all these assumptions 20 the event horizon of the spacetime 
(|3.4| ) has a simple mathematical characterization; it is simply the unique null hypersurface 
that reduces precisely to r = at late times (see || for a related discussion). 

All the metrics we study in this paper fit into the general form given by equation 
2.2 and 2.3 in |B[ with one straightforward generalization. The normalization condition 
u^Uy-q^ = — 1 of H is replaced by u^u^g^ = —1, where g^ v is the arbitrary weakly 
curved boundary metric of our paper. Following ||, it is easy to work out the precise 
equation for the event horizon within the derivative expansion. In fact it turns out that 
equation 2.16 of || continues to hold with the single proviso that all indices in that equation 
are raised by the metric g^ v rather than rf v . The expression 2.18 of ||, for the normal 
vector to the horizon, also continues to apply with the same proviso. Finally, the area form 
a on the event horizons in our paper may be defined exactly as in section 3.1 of |J. The 
formula 3.10 of || for this area form continues to apply. Lines of constant define a map 
from the boundary to the horizon of our solutions. Let s denote the pullback of a, to the 
boundary, under this map. Let (Js)^ = *s. The following slight generalization of 3.11 of 
gives a formula for the entropy current 

40SJ +1, vS™ 

(see for the definition of h). As in j|] it follows immediately from classic area increase 
theorems of event horizons in general relativity that V /i (J s ) At > 0. Consequently {J s )^ is a 
candidate entropy current with non negative local entropy production. 

We have evaluated all the abstract expressions described above for the specific fluid 
dynamical solutions presented in our paper; we present our results. In the gauge ( pT9| ) we 
find the covariantization of equation 5.4 of || 

T h = \ + \(s h + \e h + \eA (5.2) 



An explicit expression for J$ is not difficult to determine: we find 

, + Sa) 

(5.3) 



2 1 ip. otf ^ 

AGb 3 J^ = «m ( i _ 1. ( - In 2 + ^- ) -1 S 5 + — ©i* + — S b + S a 




4 

+ b' A S? 



20 It is probable that these assumptions are more restrictive than necessary.lt would be interesting to 
know which if any of these assumptions can be relaxed without affecting the simple characterization of the 
horizon that we give below. 
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where 



S a 



+ 



2 12 4 V6 6 48 3 J 12 

b 2 

— [-864 - 65 (8 + 48C + 4vr + 4vr 2 - 36(ln 2) 2 - 12vr In 2 + 16 In 2)] 

9 ^ + i+ln2 > l 61' 
4 2 



(5, 



^ T3 00 

16 B " 



1 Tjfin 

144 " 



As we have explained above, the entropy current presented in ( |5.3| ) is guaranteed by 
the area increase theorem of general relativity to have non negative divergence. This claim 
may also be verified algebraically. Using the equations of motion, we find that 

I 77 + I -hi) Ml 1 ! > 

2r] 



^ + (7r + 4 +! ln2) n^' 



+ 77 



u.c 



+ 



1 



16vrT 

, 7T 



4vrT v 2 



+ 2 + 31n 2)u a u f3 V a V (3 ^ 



(5.5) 



+ ... 



where the RHS is accurate only to third order in the derivative expansion. Of course the 
RHS, which is a sum of squares, is manifestly positive, in accordance with the second law. 

The second line in ( |5.5| ) represents the production of entropy due to the forcing effects of 
a time dependent dilaton. Several comments are in order here. Consider the forced Navier 
Stokes equations V^T^ = ft 1 . For an arbitrary function ft* this equation is consistent 
with either local entropy increase or decrease. The fact that entropy production is always 
locally non negative in our situation is a consequence of the particular form of our forcing 
function. 

In order to see what is special about this forcing function, consider a dilaton field 
that is a function only of the boundary time (see §[?] for a more detailed discussion of this 
situation). According to ( 3.20 ), the forcing function for such a field f ^ oc d^docj)) 2 . The 
sign in this expression is important; it implies that a varying dilaton field always pumps 
energy into the system. This energy raises the fluid temperature increasing its entropy. 

This qualitative feature should be true of any forcing function that results from the 
coupling of a source to a gauge theory operator. A source term for any boundary operator 
excites the expectation value of this operator away from its thermal average. The expec- 
tation value of this operator then decays over a time scale 1/T; in this process the energy 
of this excitation thermalizes, resulting in an increase in entropy. 



6. AdS 5 Kerr BlackHole 

A three parameter set of asymptotically globally AdSs rotating black hole solutions has 
been presented in |35|]. These solutions are labeled by their temperature T and their two 
angular velocities, uj\ and u>2, on S 3 . As we will see in more detail below, these solutions 
are dual to a field theory configuration which is well described by the equations of fluid 
dynamics at large T (the radius of the boundary sphere is unity in our conventions). 
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6.1 Summary of Results 

In this section we will rewrite the exact solutions of [S5| in the gauge ( 2.10p . The final 
metric admits a remarkably simple all orders re expression in fluid dynamical terms: it 



may be put in the form (4.7) with 
Wx = l ; 

W 3 = 



W 2 = l + ^ af3 u a P 



W, 



w, 



5fj,u 



2"l / 1 n R 



3r 2 a ^ 6r 2 



(6.1) 



al3 



-P, 



Moreover the radial location event horizon of this metric is given by the solution to the 
equation 



1 



(l + :2r )(r 2 H + ui)(r 2 H + u4) 



2m 



H 



The entropy of this black hole is given by 



rff (l_ w 2)(l_ w g)(l+ 1) 



(6.2) 



(6.3) 



Moreover the various components of the local entropy current defined in || evaluate, on 
this rotating black hole solution 



;2\ 



UJ7 



j4 >_ uilrjj + ui){r 2 H + 

r H {r 2 H + ul) 
"2(r 2 H + ul)(r 2 H + u 2 2 ){l - ul) 



(6.4) 



ji> 



r H {r 2 H + ujl) 



Notice that the exact result (|6.l| ) admits a remarkably simple derivative expansion. 
The derivative expansion of the functions W±, W2, and terminates at second 

order in derivatives. On the other hand the function W3 includes terms at all orders in 
this expansion: however this expansion sums up to a simple geometrical series. Several 
additional remarks are in order 



1. We check below that the expansion of (6.1) to second order in the derivative expansion 
matches perfectly with the general predictions given in (f4.7|). This is a nontrivial 
consistency check on the main result of this paper. 



Note that fO) receives contributions from terms that depend on the curvature of 
the boundary metric only at second order in the derivative expansion. Consequently 
the second order fluid results, fl6.1|) listed above, exactly account for all curvature 
dependencies in the exact result (|4.7|) . In particular the AdS Schwarzschild solution 
is reproduced exactly at second order in the fluid expansion. 
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3. The fluid dynamical expansion of the metric (O) is convergent for 2r 2 > oj a aoj a ". 
It follows that this expansion is convergent everywhere outside the horizon for T 2 3> 

Expressions for the location of the event horizon and the local entropy current may also be 
expanded rather simply, to all orders, in fluid dynamical expansion. We demonstrate below 
that the expansion of these expressions to second order also agree with the predictions of 
the previous section. 

In summary the fluid expansion works even better for the exact solution ( |4.7| ) than 
we might naively have the right to expect. In the rest of this section we will derive and 
explain all these results in detail. 

6.2 Rotating Black Holes in Fluid Dynamical Coordinates 
6.2.1 The AdS Kerr Solution and its Stress tensor 

The asymptotically globally AdSs rotating black hole solution presented in (6!| is given by 



A r / ~ ui sin 2 9 2 u 2 cos 2 6 A p 2 2 p 2 

—n I dt 2~# —dip) + —dr +— - 

P 2 V 1 — ujf 1 — ^2 / 



2 , P jcv2 



A / / r 2 +uj 2 \ 2 / r 2 -\-uj 2 \ 2 \ 

+ -J f S in 2 e ^di - - r^ #j + cos2 y^dt - - ) ] (6.5) 



1 + r 2 / / r 2 +oj 2 A / r 2 + w 2 - u2 
H ^-2" ( u>2 s ^ n2 ® ( \d(\> J + wi cos 2 I u>2di ^-d^ 



where 



2 2 , 2 2 r\ , 2 ■ 2 r\ 

p = r + u>i cos 9 + u>2 sin 6 

A r = 1(1 + r 2 )(r 2 + o, 2 )(r 2 + l^ 2 ) - 2m (6.6) 
Ae = 1 — uj\ cos 2 — w| sin 2 

In the large r limit (i.e. near the boundary) the induced metric on a surface of constant 
r is given by 



r/s 2 

Bnd " g^dx^dx" 



A e di 2 dQ 2 sin 2 © ~ jf . 2 cos 2 . , 7 lfs2 

(1 — tt»f)(l — ) Ae 1—Ujf \—i02 



(6.7) 



Consequently the solution ( p.5|) is dual to a state of the CFT on the space Weyl equivalent 
to ( |6.7D . (6.7) describes the so called Rotating Einstein Universe. 
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The boundary stress tensor ( 2,13 ) dual to this metric is easily evaluated, and may be 
written in the form 



V = g^(<V + ^u v ) + ^JL- (r^R^ _ ^ ) where 



-Undx^ = di — Ul S ™ ® rf<A — ^2 cos Q ^j; 

1 — 1 — or 



(6-8) 

A^di u\ sin 2 * - ^2 cos 2 



, .,. 2 -(d^ + widt) —(dip + u 2 dt) 

(l-wf)(l-w|) 



" <9i 

where R^ v \a denotes the boundary curvature tensor associated with the metric g^. Apart 
from an additive curvature dependent piece whose form is dictated by the conformal 
anomaly 21 , this stress tensor is precisely that for a perfect conformal fluid with velocity 
vector u M = (1, 0, 0, 0) and pressure p = %™ Gb ■ 

6.2.2 Relationship to analysis in ^ 

We pause here to connect these results to those derived in ||. The authors of that paper 
analyzed the same rotating black hole solutions, however they worked with a different set 
of coordinates denoted here by (f,t,(f)',ip',@'). In these coordinates the boundary of the 
rotating black hole solutions is naturally S 5 x R as we now explain. 

The Rotating Einstein Universe is Weyl Equivalent to S 3 x R. In order to see this 
consider the boundary coordinate transformation 

2 „ (1 - ul) cos 2 0' 

COS Z = \ ^ 

1 - to 2 sin 2 0' - loI cos 2 0' 
. 2 A _ (l-a; 2 )sin 2 0' 

Sm l-u; 2 sin 2 0'-u; 2 cos 2 0' ( 6 ' 9 ) 

<j> + uj\t = 4> 



Expressed in terms of these new variables, the boundary metric ( |6,7| ) may be written as 

^ = ^-T7^ (~ dp + d& ' 2 + cos2 '^ /2 + sm 2 0'#' 2 ) 

H l-ujf sur 0' - oot, cos 2 0' V / 

= 7 2 (-dP + dn§) ( 6 .10) 
1 



'l-u; 2 sin 2 0' - uj cos 2 0' 
where d£l\ represents the usual metric on the round 3 sphere. 



21 Note, for example, that the trace of the stress tensor in this background as calculated from above |67j 
is 647^5^ = ^/3 R al3 -Si. 
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In order to obtain a boundary metric that is actually that on S 3 x R rather than 
simply conformal to it, the authors of p| worked with the redefined radial variable r = 

, J- = , n = n so that the induced metric on slices of large 

y/i-Lj2sm 2 e-w% cos 2 e' cos2 ©- w l sin © 

constant r is equal to 

dn 2 

?± r = -dt 2 + dni (6.11) 

With these conventions the rotating black hole solution is dual to a fluid flow on S 3 x 

time. The fluid velocities and temperatures of this flow are given by acting upon the velocity 

i 

and temperature fields = (1,0,0,0) T = QMH ; Q f the previous subsection, with the 
coordinate transformation ( |6.9| ) followed by a Weyl transformation (compare ( |6.10| ) and 
(|6. 1 1[) ) yielding 



( d d d \ 

f = ^A 7 (6.12) 

7T 

1 

7 = 



'l-ul cos 2 8-a;f sin 2 9' 
precisely as reported in 0. 

6.2.3 Recasting the solution in the Fluid mechanical gauge 



Continuing with our analysis of ( |6.7| ), we will now recast this metric in the gauge ( 2. ICQ . 



As was explained in || , the coordinates in gauge ( [2.10D are adapted to a congruence of null 
ingoing geodesies, whose tangent vectors near the boundary are given by 4| = 1, = 
in Graham Fefferman coordinates. While it turns out to be difficult to solve for the most 
general null geodesic in (|T7]), precisely this congruence turns out to be easy to determine 
and is given by 



dr 
dX 


= -1 




de 

~d\ 


= 




dt 


(r 2 


+ u\){r 2 + ojl) 


dr 




r 2 A r 


d4> 
dr 


(1 

= ^ 2 


- u 2 ) dt 
+ ujf dr 


di> 
dr 


(1 

= U 2 ~^ 


- oj 2 ) dt 
+ uj 2 dr 



(6.13) 



As we have explained before, in the gauge ( 2.1C| ) the coordinates are constant along 



this congruence of geodesies, and the coordinate r is simply the affine parameter along 
these geodesies. Consequently, it follows that the change of variables to the Eddington- 
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Finkelstein like co-ordinates ( [2.10 ) is given by 



dt = dt- - — V — —^-dr 

r 2 A r 

d--?)(r 2 + -i), r (614) 



r 2 A r 
{r 2 +UJ 2 )(1 _ UJ 2 ) 

r 2 A r 



Expressed in terms of the new coordinates the metric becomes 

t 2 _ { - u>x sin 2 G , u>2 cos 2 8 , A 
ds = 2dr dt ^— do s— dip 

-%(dt- " lSm2 2 % - " 2CQS2 2 9 #1 2 + ^0 2 

p 2 V l-^i 1-^2 / A © 



(6.15) 



A / / r 2 +uj 2 \ 2 ( t 2 + lo 2 N 2 ' 
H „- sin 2 6 I u}\dt ^-d^ I + cos 2 I uj2dt §-dip 

P \ V l-UJ 1 J V l~ w 2 

j_ -|- / / r 2 + lo 2 \ ( r 2 + a; 2 
H 2—2- I W2 sin 2 I widt ^d(p I + tc>i cos 2 I uj2dt ^-dip 

As we describe in detail in the appendix |B|, we conclude that the AdSs-Kerr metric 
can be written in the following manifestly Weyl-covariant form 

ds 2 = -2u^dx^(dr + rA\dx x ) + ^g^dx^dx" 

1Z \ 2m ( 1 \ — 1 

U[jV\uj X v + uj^ujxu + —u^Uu j dx^dx u + — ^ ( 1 + J u [1 u v dx >i dx v 

= —2WiUf M dx >1 (dr + rA\dx x ) 

+ [r 2 (W 2 g^u + Wzu^Uv) + r (W^u v + u^W iu ) + W 5 ^] dx^dx v . 

(6.16) 

with the functions in these equations given by (|6.1|). 
6.3 Horizon and Entropy Current 

As we have mentioned above, the location of the event horizon (in the radial variable of the 
gauge ( [2,10D ) is given by the solution to the equation (|6.2|) , The solution to this equation 
is easily obtained at second order in the derivative expansion: we find 

1 ( b \ 2 2- 

- = (2m) 4 
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It is easy to check that this result agrees with our general prediction ( p,2| ), once that 



equation is re-expressed in the radial variable of the (2.10) 



1 h f 1 1 \ 37> 2 r°° 

r *=s + i(*n 6B n 6i 'WjL* <Jr) ( " 8) 

On the velocity configuration dual to the black hole 65 and 65 vanish and 

UJ a /3iO /v / U a gU 

Hir) — 



6 ./„ w 4r 



inputting these values we see that (|6l|) and (|07|) agree 



In a similar fashion, the formula for the entropy and entropy current ( |6.3D for these 
black holes may be expanded to second order: we find 

4 \ 2 6 



(6.19) 



s = 


n 3 


& 3 (i-^ 2 )(i- 


J l = 




J* = 




b 




£J 2 (1 - UJ 2 2 ) 


b 



On the other hand the general formula for the entropy current ( |5.3| ) evaluates, on the 
specific fluid flow at hand, to 

4G 5 6 3 Jf = u» (l-j (5u a pu a P + ^P^V x lu„ x (6.20) 



It is easy to verify, using ( |B.4[) , that ( |6.20D and ( 6.19j ) agree. Upon integration S = J i/gJ 
with y/g = we also reproduce the first of (gig) . 

7. Some Simple Solutions of Forced Fluid Dynamics 

In this section we construct some simple solutions to the equations of motion of fluid 
mechanics with forcing terms derived above, and consider their interpretation in the bulk. 
In order to make the physical points of interest most immediately, in this section, we often 
work with the crudest approximations that capture the physics at hand. We postpone a 
more careful study of these solutions (and hopefully several others) to future work. 



In subsections 7.1 we study fluid flows that are pushed to high Reynolds numbers by 
the effective forcing effect of a varying background metric. In subsection |7.2j we study 
'cosmological' solutions of fluid dynamics corresponding to a time dependent but spatially 
homogeneous dilaton. As we have explained in the introduction, we believe that these 
solutions qualitatively capture the excursion of the bulk geometry into regions of strong 
curvature. Finally, in subsection ffl^ we study a fluid that is pushed into motion by the 
forcing effect of a varying dilaton field. 
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7.1 Metric as the Forcing Term 

7.1.1 Hydrostatic Solution in an arbitrary spacetime 

We begin by considering the case where the forcing term arises due to a non-trivial metric. 
We start with the stress tensor in the perfect fluid approximation, 

= {nTfig^ + AuV), (7.1) 

in the presence of a time independent metric of the form, 

ds 2 = goodt 2 +g ij dx i dx j . (7.2) 

It is easy to see that the conservation equations, 

V M T^ = 0, (7.3) 

then admit a hydrostatic solution, with four- velocity, 



u ° = V^oo, ij = 0, (7.4) 

and temperature, 

V l^oo I 

The temperature dependence can be understood as arising due to the gravitational red- 
shift. We see that the temperature is higher, where |goo| is smaller, i.e., where the grav- 
itational potential is deeper. In the bulk the radial location of the horizon is given by, 
r = 7rT, and becomes a function of the spatial coordinates. 

This solution gets corrections due to the additional terms in the stress tensor, eq.( 2.13| ). 



Since the fluid is at rest the viscosity term is irrelevant. Although we have not attempted 
this, it should be possible to work out the corrections due to the second-order terms found 
in this paper in a straightforward manner. All such terms are suppressed by a factor of 
(Ti?)~ 2 , where R is the length scale associated with the boundary spacetime, compared 
to the leading order result presented above. These terms will lead to a corresponding 
second-order correction to the dual geometry. 

7.1.2 Small deviations from hydrostatic equilibrium 

We now specialize to a boundary metric of the form 

5oo = 5oo (z) 

{z is one of the spatial coordinates) and study a particular perturbation about hydrostatic 
equilibrium. In particular, we search for a steady state solution in which, in addition to 
a varying temperature, the fluid has a small, time independent velocity purely in the z 
direction. Throughout this subsubsection we work to first order in this spatial velocity u z . 
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The v = component of eq. ([7.3|) gives 22 

rV|ffoo| = c, (7.7) 



where c is an integration constant. The v = z component of eq.( [7.3|) is unchanged from 
the hydrostatic case, at linear order, and gives, again, eq. (|7.5|) . From, eq. (|7.7|) , (|7.5|) it 
follows that 

u z =u z o\goo\, (7.8) 



where Uq is a constant. ( [7.5|) , and ( |7.8| ), determine the temperature and velocity T,u z , 
as a function of z. Let us suppose that goo = — 1 at a particular location and changes 
to — (1 + /) over the length scale L. It follows that the change in the fluid velocity over 
the same length scale equal in magnitude to fuff. While Uq has been assumed to be small 
compared to unity in this subsection, and while we might wish to restrict / also to be 
small compared to unity for some physical purposes, it is consistent to hold each of these 
quantities fixed as TL is taken large. In this case the formal 'Reynolds number' of this flow 

R e ~ AvTL, (7.9) 

becomes parametrically large in the fluid dynamical limit TL S> 1. 

While the static flow described above is rather tame, we believe it illustrates the general 
point argued for in the introduction, namely that even though the forcing effects of gravity 
are mild (suppressed by 1/L) they can build up to 0(1) changes in the velocity - and hence 
high Reynolds numbers - over length scale L. Consequently one might hope to be able 
to stir the fluid into steady state within the validity regime of the approximations of this 
paper, though we would probably need a more general metric (one that depends on several 
spatial directions as well as time) for this purpose. 

7.1.3 Flows forced by a time dependent spatial metric 

We can also consider an analogous situation where the metric is, 

ds 2 = -dt 2 + g zz {t)dz 2 + dx 2 + dy 2 , (7.10) 

with only one non-trivial spatial component that depends on time. A consistent solution 
can be obtained with u°,u z being the only non-zero components of the 4- velocity. Working 
with the perfect fluid stress tensor, one finds from the v = component of eq. (|7.3| ) that, 

rVV = t\f2- ( 7 - n ) 

{9zzY 12 

The four- velocity satisfies the constraint, 

(u ) 2 - {u z fg zz = 1, (7.12) 



22 We work in the perfect fluid approximation and at first order in velocities, ft should be straightforward 
to account for the first nonzero corrections to these approximations. 
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The last equation is the v = z component of eg. (7.3). Using eq. ([7.11 ), eq. fl7.r2D , to elimi 



nate two of the variables, one can reduce this problem to quadrature. We skip some of the 
details here. 

The temperature and 4-velocity are time dependent in this case. It is easy to see that 
again order unity changes in the metric, g zz , result in order unity changes in the velocity. 
The time varying metric allows energy to be pumped in or out of the system. If the time 
scale over which this happens is big compared to the temperature, as is needed for the fluid 
mechanics approximation to be valid, one finds analogous to the the spatially varying case 
discussed above, that the Reynolds number is much bigger than unity. 

7.2 Time Dependent Dilaton and Highly Curved Spacetimes 

We now study a simple 'cosmological' solution of fluid dynamics. In the solution we study 
the metric is taken simply to be tw, but the dilaton is taken to be a slowly varying function 
of time. We are interested in answering the following question: suppose we start with field 
theory heated up to a temperature T and at coupling constant Ai. Suppose we then slowly 
vary the coupling from Ai to A2 over a time At. What is the final state of the theory at 
the end of this process? 

Using the equations of fluid dynamics derived in this paper, it is easy to answer this 
question (we comment below on the validity of these equations for the purposes of this 
issue). In particular, there exists a solution to the fluid equations of motion in which the 
fluid is always at rest and the temperature remains spatially homogeneous but heats up 
slowly in time, in response to the varying dilaton. In more detail consider the configuration 

u° = 1, u l = 0. (7.13) 

Taking the leading contribution to the dilaton forcing function (this is the first term on 
RHS of the second equation in eq.( [2.13D ) gives, from eq.( |3,20| ). 

V M T^ = -(7rT) 3 <9>d o 0- (7.14) 



If we take the perfect fluid stress tensor, eq.( [7.1[) , and 4-velocity, eq.( 7.13j ), this gives rise 
to one non-trivial equation for T, which can be integrated to give, 

T(t f )-T( tl ) = ^- f\j>?dt. (7.15) 

We see that, irrespective of the details, time dependence of the dilaton always increases 
the temperature and thus the energy and the entropy density of the fluid monotonically. 
In the bulk correspondingly the horizon area increases monotonically. This ties in with the 
discussion below (|5.5| ) . 

Note that 

T TAt v ; 

i.e. a total change in the dilaton, Acfr, over a time, At, results in a small fractional change 
in temperature if 5t is sufficiently large. When the LHS of ( |7.16| ) is small the condition for 
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the validity of fluid dynamics 

f (A0) 2 



< 1. (7.17) 



T 2 (TAt) 2 

is automatically met provided A<f> > 1 and hence TAt > 1. 

In summary a time dependent coupling constant always heats up the system; however 
the fractional increase in the temperature can be arranged to be small - even for a large 
fractional increase in the coupling - provided the coupling constant is changed relatively 
slowly. 

Corrections to eq.( 7.15| ) will arise because of corrections in the stress tensor, eq.( 2.13j ), 



and in the forcing function, eq.( |3.20| ). If the condition, eq.( 7.17| ), is valid corrections due 
to second order terms in the stress tensor are small 23 . Corrections in the forcing function 
go like, T 2 cj)(j). These are small, compared to the leading order term, if TAt > 1. 

7.2.1 Highly Curved Spacetimes 

In the bulk of this paper we have derived a set of fluid dynamical equations from the 
Einstein Hilbert Lagrangian and studied some aspects of the dynamics of these equations. 
Our analysis so far has been quantitative. On the other hand, in this subsubsection we 
attempt to explore some qualitative aspects of bulk dynamics of highly curved spacetimes 
via fluid dynamics. 

In particular, we wish to investigate the behavior of a set of classical solutions of string 
theory (rather than gravity) that is dual to evolutions of large N c N = 4 Yang Mills 
theory with a time varying 'tHooft coupling that is lowered down to order unity. The 
bulk equations that govern these systems are unknown. Nonetheless we will attempt to 
qualitatively understand these evolutions using an expectation based on physical intuition 
: at long wavelengths these (unknown) equations of classical string theory should reduce 
to the equations of boundary fluid dynamics with A dependent parameters. Though we do 
not know the detailed A dependence of these parameters, we will attempt to estimate their 
qualitative properties to the extent needed for our analysis. 

We pause here to explain this intuition in more detail. The equations of fluid dynamics 
follow simply from symmetries combined with the physical expectation of local equilibra- 
tion. Any interacting system is expected to equilibrate locally over a length and time scale 
set by its 'mean free path' which may be crudely estimated by the ratio r//p where r\ is 
the viscosity and p the energy density of the system. It follows from 'tHooft scaling and 
dimensional analysis that this ratio takes the form q(X)/T for some function q(\) in AA=4 
Yang Mills theory. We know from perturbation theory that q(X) diverges at weak coupling, 
and we know from the AdS/CFT correspondence that q(\) is a constant at strong coupling. 
From these two behaviors, it seems reasonable to guess that q(X) decreases monotonically 
as the coupling is increased or, at least, that it never diverges at a finite value of A. If this 
last guess is correct, the effective equilibration length scale for M = 4 Yang Mills should 
be of order 1/T for all A greater than or of order unity, and fluid dynamics should be a 
good description of any evolution of the system in which all quantities vary on a length 



3 Any corrections due to the viscosity term vanish in this case. 
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scale that is large compared to 1/T.We make this assumption in what follows below. We 
also assume the absence of any phase transition in high temperature N = 4 Yang Mills as 
A is varied from large values to values of order unity. We wish to emphasize that if either 
of the assumptions listed in this paragraph fails, none of the conclusions we reach in this 
section need apply. 

With all these caveats in mind, we now proceed to consider a situation where we start, 
in the far past, with a value of the dilaton on the boundary such that gsN ^> 1. We take 
the dilaton to be a slowly varying function of time and reduce its value till, g s N = \ m in- 
Then increase it again so that in the far future g s N S> 1. 

The fluid dynamical solution presented above accurately captures the boundary dy- 
namics of this situation provided A min 3> 1. When this condition is met, the results of this 
paper eq.( |3.4[ ) also yield the bulk dual of this field theoretic evolution. 

Now let us instead consider take X m i n to be 0(1). In this case we expect the evolution 
described in this subsection to be given by some classical solution of string theory rather 
than supergravity. We also expect spacetime curvatures to become of order string scale - 
and so for stringy effects to become very important - over the times at which A is of order 
Xmin- However curvatures are small in string units (and so the SUGRA approximation is 
good) at early and late times. In other words the dual bulk solution is expected to be a 
cosmology where an initially smooth space-time becomes highly curved, with a curvature 
of order the string scale so that it is not well described by the two derivative SUGRA 
approximation, and then returns in the future to being smooth again. 24 We now want to 
inquire what information we can glean about this evolution process by the use of boundary 
fluid dynamics. 

Of course the specific fluid dynamical equations we have derived in this paper are only 
valid at very large A. However aspects of these equations are dictated merely by symmetries 
and the physical expectations of equilibration, and may have much greater validity. Of 
special importance to us in this section is the structure of dilaton dependence in the forcing 
term in the Navier Stokes V pT^ v = f u studied in this paper. As we have explained in the 
introduction, it follows from general field theoretic reasoning that = e'^Cd^cf) where 
C is the expectation value of the Lagrangian. As we have explained above, it further 
follows Weyl invariance that upto the one derivative level e~^C = T 4 r(A) + T^g(X)u.d4>. 
Similar statements apply to the form of the stress tensor upto first order in the derivative 
expansion. It follows that the leading two (boundary) derivative equation of motion that 
governs the evolution described in this subsection is 



where h{\) is an unknown function that evaluates to unity at A = oo and is expected to 
be of order unity provided for A greater than or of order unity. On the other hand g{\) is 
another unknown function such that g(X) — » as A — > oo. There seems no reason not to 
assume that g(\) is of order unity when A is of order unity. 

24 For some recent work on time dependent cosmologies in the context of AdS/CFT, see ItS, Ffj, [77], F3, 




(7.18) 



§, §, §, 
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Consider again varying the coupling constant down from infinity. The first term in 
fl7.18|) dominates the evolution of the temperature when the coupling is large enough (i.e for 
g(X) <C T~ 1 do4> ) an d ( 7.15| ) applies in this domain. On the other hand once g(X) increase 



above T 1 do4>, the second term in (7.18|) dominates the evolution of the temperature. The 



details of the subsequent evolution depend on the precise form of the unknown function 
9(A). 

In order to get a sense for how this works, however let us assume that g(X) = K/X a for 
some positive power a (this would be the expectation for g(X) at large A if it was generated 



at finite order in the a' expansion). In this case the second term of (7.18) dominates over 
the first when A < Ao ~ (^) 1 ^ a ' ^ or lower values of A we can ignore the first term in 
(7.18), and then that equation is easily integrated. We find that, in lowering A from Aq to 



X m i n , the temperature changes by a factor e a x min A o . Now ^ must be small in order 
for fluid dynamics to apply. It follows that Ao 3> A m j n when \ m in is of unit order, so that 
lowering A from Ao to A m j n changes the temperature of the fluid by a finite factor that is 
is independent of Ao- 

In summary, if the coupling is lowered very slowly from A to Ao the temperature 
is almost unchanged. Subsequent lowering of the coupling from Ao to \ m i n changes the 
temperature by a fixed finite factor that depends on X m i n but not on Ao or other details 
of the process. We conclude that the principal conclusions of the previous subsection - 
namely that the dilaton may slowly be lowered from an arbitrarily large though fixed value 
to a value of order one with a finite (though, in this case, not arbitrarily small) change 
temperature - carry through unchanged. This appears to be an interesting statement about 
the dual bulk evolution. We have argued that nothing particularly dramatic happens to a 
spacetime - at least if it is at large finite temperature - as it traverses through a region of 
string scale curvature. We emphasize that it was very important to our analysis that we 
were at high temperature; the fluid dynamical analysis presented here sheds no light on 
the evolution of the same bulk geometries at zero temperature. 

In fact one can try and go even further. After having reduced the dilaton to a value 
where, g s N ~ 0(1), we can reduce it even more till g s N <C 1. At this stage the notion of 
spacetime geometry has completely broken down and the correct description is in terms of 
a perturbative Yang Mills theory. Thereafter one can increase the dilaton, till in the far 
future it meets the condition, g s N 3> 1. As long as the driving force due to the dilaton is 
slowly enough varying, one expects that the boundary theory will be well defined, and the 
initial state we start with can be continued past the highly curved region to match with a 
smooth geometry in the future. 

Now fluid mechanics can continue to be a good description, at small g s N, provided, 
the function h(X) does not blow up 25 as A — > 0, and, the temperature varies slowly enough, 
satisfying the condition, 

Imfp^ < 1, (7.19) 

(where l m f p stands for the mean free path). The condition ( 7.19 ) is expected to translate 



3 This seems very reasonable and can be checked against a perturbative calculation. 



- 39 - 



into 

^2 « A (7.20) 
at small A, a stronger condition than ( 7.17] ) which holds at A of unit order or greater. (This 



expectation follows from the estimate, that in perturbative Yang Mills theory 

lmfp ~ J^NfT- (7 - 21) 

as opposed to l m fp ~ y at strong coupling) . Consequently, requiring that fluid mechanics 
is a good effective theory in the perturbative regime puts more stringent limits on how fast 
the dilaton can be varied. Prom eq. (|7.19| ), eq. (|7.21| ), we get, 



< T. (7.22) 



dt 

Nonetheless, using a dilaton profile which varies slowly enough to meet this more stringent 
condition, we can construct solutions where two asymptotically smooth regions of spacetime 
are connected by an intermediate region which is highly curved (dual to Yang Mills theory at 
any specified - though nonzero - coupling) . This intermediate region admits no conventional 
description in terms of spacetime, but we expect it to be well described by fluid mechanics 

26 

As must be clear to the reader, the discussion this subsection makes several assumptions 
and is tentative in some respects. In particular we reiterate that it makes assumptions 
that are difficult to directly verify: in particular the assumption that there are no phase 
transitions in high temperature N = 4 theory as g s N is varied from large values to 0(1), 
or smaller values. 

7.3 Spatially varying Dilaton 

Finally we consider a spatially varying dilaton as a forcing function. For simplicity we 
take the dilaton to depend on only one spatial coordinate, x. We also set, the u y ,u z , 
components of the four-velocity to vanish. As in our discussion of the time dependent case 



above, we start with the perfect-fluid stress tensor, eq.(7.1), and leading dilaton forcing 



function. Eq.( [7,14| ) then gives rise to two condition. From, the v = component we get, 

d{T A u°u x ) 



dx 

And, from the v = x component, 



0. (7.23) 



rf (TW)3 + 1) = _r3 

dx 7T 

where prime indicates derivative with respect to x. In addition the condition, 

(-u ) 2 - {u x f = 1, (7.25) 



26 Also it is worth noting that once the dilaton is small enough, we could more directly analyze the time 
evolution in perturbation theory, even when fluid mechanics is not valid. This may even allow us to study 
evolutions that continue to strictly zero field theory coupling in finite time. 
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is also valid. This yields three equations in the three variables, T,u°,u x . Eliminating two 
variables from eq.( 7.23j ), eq.( 7.25 ), and substituting in eq,( [7.24 ) reduces the problem to 
quadrature. 

Here let us consider the case where where the spatial velocity is small, i.e., 



u x < 1, 



(7.26) 



and work, to begin with, to linear order in u x . To this order, u° = 1, and from eq.(7.23), 

u x = (7.27) 

where c is a constant. Eq,(|7gg), now gives, 

T 5 (x) 



i_ JLl£ 



k'\2 



ydx. 



(7.28) 



where Ti,uf stand for the temperature and velocity at x — » — oo. Without loss of generality 
we can choose conventions so that uf > 0. Then we see that as x increases, the temperature 
decreases and the velocity increases. From, eq.( 7.27 ), eq.( 7.2Sj ), we see that the length scale 
of variation of the temperature and velocity, are set by the forcing function of the dilaton. 
We denote this length scale by L below. 

The analysis can be improved by working to quadratic order in u x . In general also the 
dilaton will depend on all spatial coordinates, this will lead to spatially varying velocities 
and to viscosity effects being important. The Reynolds number for a fluid flow is given by, 



Re 



Au x TL, 



(7.29) 



where Au x is the variation in the velocity along the flow. From, eq.(7.27), eq.(7.28), we 
see that, for A(f> ~ O(l), Au x ~ wf^f- ~ (uf) 2 jf^. This gives Reynolds numbers of order 
unity or smaller. Again this is in line with the expectations spelt out in the introduction. 
The smallness of the Reynolds number of this flow is connected to the fact that the dilaton 
yields a forcing function that is of second rather than first order in derivatives. 

Corrections to these flows can be systematically computed, using the corrections to 
the stress tensor, eq. fl2.13 ), and forcing function, eq.( 3.20| ). 



8. Discussion 



In this paper we have used the AdS/CFT correspondence to determine the metric dual 
to an arbitrary flow of the fluid of a conformal field theory on an arbitrary weakly curved 
four manifold, with an arbitrary slowly varying coupling constant. We have also explicitly 
constructed the event horizon of these metrics in the derivative expansion, and used grav- 
itational ideas to propose the construction of a dual entropy current whose local increase 
is guaranteed by the area increase theorem for black holes in general relativity. 

We have tested the constructions of this paper against an exact class of asymptotically 
globally AdSs solutions of Einstein's equations, namely rotating black holes in global AdSs. 
We have demonstrated that these solutions admit a remarkably simple rewriting in fluid 
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dynamical terms, and have verified that they agree in detail with the general constructions 
of our paper when expanded out to second order in a derivative expansion. 

The connection of large rotating black holes to fluid dynamics was previously partially 
explored in [g]. The authors of that paper noted that the stress tensor dual to rotating black 
hole solutions was exactly reproduced, to appropriate order in the derivative expansion, by 
the predictions of fluid dynamics. However these authors also noted an apparent first order 
discrepancy between the properties of charged rotating black holesp7], p^, 7C, 89] and the 



predictions of charge fluid dynamics. A generalization of the calculations of this paper to 
account for fluid charge would allow us to pinpoint the source of this worrying discrepancy. 

Turning to another issue, we believe that it should be possible, within the approxima- 
tions of this paper, to choose a time dependent fluctuations of a flat space metric to stir 
the conformal fluids studied of this paper into steady motions with high Reynolds num- 
bers. Even using the forcing functions provided by linearized fluctuations of the metric 
away from flat space (see the introduction) we believe it should be possible, for instance, 
to stir a fluid on R 3 or the sphere into configurations that resemble classic experiments 
that display turbulence 27 . Consequently, the map from fluid dynamics to metrics in this 
paper should yield the spacetime dual to a turbulent fluid flow. The understanding of ap- 
parently universal exponents in turbulent motions remains one of the outstanding unsolved 
problems of statistical physics. It is conceivable that the dual spacetime perspective could 
permit qualitatively new insights for such flows. 

Of course transient turbulence occurs in fluid flows even without forcing. For instance, 
a configuration of our conformal fluid at temperature T on an S 3 of radius R, and with 
(varying) velocities of unit order has Reynolds number of order TR and so presumably 
undergoes turbulent motion. Now all such fluid configurations eventually settle down into 
(presumably unique) non dissipative solutions of of fluid dynamics with given conserved 
energies and angular momenta; the rigid rotations studied in Q and in section ||| of this 
paper. Let us translate this expectation to bulk language. Consider an arbitrary gravita- 
tional configuration, in global AdSs, with appropriately large energy and angular momenta. 
Such a configuration will undergo gravitational collapse and eventually settle down into 
one of the rotating black holes studied in ||]. The gravitational dynamics of this settling 
down process could be dual to a turbulent fluid flow. It may be interesting to pursue this 
connection further. 

Finally, we end this paper with an amusing thought. Recall that we have been able 
to construct a bulk dual to fluid flow in an arbitrary boundary metric. Now the boundary 
metric can, in particular, be chosen to have a horizon; for example it could be taken to be 
the metric of a 4 dimensional Schwarzschild black hole(See |)0| for an example). It may be 
interesting to investigate the bulk dual description of the boundary fluid falling through 
horizon of this metric. 
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7 A time dependent hot provides a forcing function in the i th direction. 
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A. Force on the Boundary fluid due to a Varying Dilaton 

In this appendix, we calculate the force that a varying dilaton field applies on a fluid at 
the boundary. To that end, we start with the bulk expression for the boundary energy- 
momentum tensor 



16vrG 5 T^ = lim r 4 {2{K CD h CD h^ - K 



r— >oo 



B) 



1 / h A h CD \\ ( A -! 

+Qb ~^ A b- \ (V A $V B $ - VcW V * 

where HaBiKab and Gab are respectively the induced metric, the extrinsic curvature 
and the Einstein tensor of the constant r hypersurface. V is the covariant derivative 
corresponding to the constant r hypersurface in the bulk. 

To compute the divergence Jb = V^T^ 4 , we proceed as follows - from Gauss-Codazzi- 



Mainardi relations (See, for example, Eqn. 10.2.24 of |91]), we have 



2V A {K CD h CD hi - K%) = -2R CD g%n D = -V c $V D %gn D = -V B $V n $ (A.2) 

where in the last step we have used the Einstein equations in the bulk. Further, we also 
have 

1- /- - h A h CD - \ 1- - 

--Va V A W b $ - Vc$V D $J = _-V 2 $V B $ (A.3) 

Using these two equations along with the reduced Bianchi identity V aQ AB = 0, we 

get 

167rG 5 V A Ti = -V B (v n <5> + ^V 2 ^ (A.4) 

Multiplying by r 4 and using IQttG^C = - lim^oo r 4 (V n $ + \V 2 <&) we get the follow- 
ing relation among the boundary variables 

V M T£ = e-*CV v <t> (A.5) 
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B. AdS Kerr 



In order to rewrite AdS Kerr metric in the fluid dynamical form, we find it convenient to 
re-express it as 

*» = 2dT U _ _ SS£!^*) - p + 1) A °f 

+ /) 2 ^- + (r 2 + uf) f n ® (# + w l( it) 2 + (r 2 + ^l ^os 6 + ^ 2 
Ae 1 — uf 1 — 

2m / sin 2 G u;2 cos 2 © \ 2 

( uj\ sin 2 G u>2 cos 2 \ 

= 2dr at ^-dcp 5— dtp 

\ 1 — u>{ 1 — ^2 / 

A 6 (ft 2 dQ 2 sin 2 8 /I( ,, 2 cos 2 G /I( , ,,\ 

71 2T71 2T + a- + 1 2# + ^ + 1 2 W + ^ 

(1 — )(1 — Wg) Ae 1 — ujf 1 — ^2 / 

^@(it 2 , o 2 <"> 2 ■ 2^\dG 2 

9T — h (w-i cos 8 + WnSm Gj — — 

.. -J 2 )(l-^ 2 2 ) A e 

a>? sin 2 6 , , , , n9 w 2 cos 2 G,,, 

+ — - n- (# + CJldi)^ + — n- (# + I 

1 — Ujf 1 — (x>2 



r 2 - 



A @ dt 2 



1 2 

2m / wi sin 2 G cl>2 cos 2 Q 
+ T 2 " ^ " " T^F # ; 



(B.l) 



This metric may be rewritten as 

2m 



ds = —lu^dx^dr + r g^dx^dx" H ^u^u u dx^dx u + Ti^ v dx^dx u (B.2) 



where 



3 = - A ^ t2 + ( w 2 cos 2 G + a, 2 sin 2 Q)^ 



w 2 sin 2 9.,, , ,n w 2 cos 2 G /I( 
+ 5-(d0 + widt) 2 + -f =-(d^ + t^dtV 

1 — Ldf 1—^2 



(B.3) 



The first three terms in flB.2j ) are simply the ansatz (2/7) while the last term, proportional 
to represents the derivative corrections (|2.8| ) to that ansatz. 

It is possible to express £ MJ , entirely in terms of Weyl-covariant curvatures and Weyl 
covariant derivatives of the velocity 
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The non-zero Weyl covariants associated with this fluid configuration are 

Afi = ; Opv = ; C^xa = 
TZ = 6 (l + uj\ + uj\ - 3(u; 2 cos 2 6 + uj\ sin 2 6)) 
uj^uj^ = 2{uj\ cos 2 + uj\ sin 2 9) 

-uJu^dx^ A dx v = sin cos OdO A ( — UJl 9 d<p ^ 2 9 dip) 

2 \1 — ujf 1 — u>2 J 

V X UJ X ^dx^ = P^V X UJ X u dx^ + UJ^UJ^U^ 

= -2dt{uj{ cos 2 O + ujj sin 2 9) 

/ uj\ sin 2 , , ujo cos 2 , , \ 

+ 2H 2 - a; 2 ) sin 2 cos 2 f ^ - 

\1-uj 2 1-ujV 

i , ,, ( wi sin 2 , , aj2cos 2 0,,\ 
P»V x uj x u dxV = -2 -J „-# + -f 

+ 2(^ 2 - u, 2 ) sin 2 cos 2 ( ^ - 

n , u i v d<d 2 sin 2 2 cos 2 , , , 

Paudx^dx" = —— + „-# 2 + 

Ae 1 — ujf 1 — ^2 

. 2r, 2m/ w 2<# \ 2 

— sin cos 7T ^ 

\l-uj 2 l-uj 2 ) 
dQ 2 

uJ l j X uJ\vdx^dx v = —{uj\ cos 2 + uj\ sin 2 0)-r — 

— Ae sin COS t: o 

V 1 — uji 1 - ujZ I 



(uj, X uj Xv + ^^P^) d^dx" = -\{ujI cos 2 + uj 2 sin 2 0) ^ 



A© 



(B.4) 



2. 9 9 . 9 _ / sin 2 , ,9 cos 2 , , 9 N 

+ -(W 2 COS 2 + ^2 Sin 2 0) 2# 2^ 

3 y l — uj^ i — UJ2 j 

A e + 2 . 2 ^ 2 n / W 2# ^ " 

H sin cos 1 



1 — u/-^ 1 — UJ2 

Using the above expressions, one can rewrite £^ completely in terms of the Weyl- 
covariant derivatives of the fluid velocity and Weyl covariant boundary curvatures: we 
find 



T lpLV dx^dx u = - ^u^T> x uj x u + uj^ujxu + ^u^uj] dx^dx v 

and p 2 = r 2 + ^-uj^uj^ = r 2 (l + ^uj a puj a ^j 



(B.5) 
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C. Notation 



We work in the mostly positive, ( — h + . . .), signature. The dimensions of the spacetime 
in which the conformal fluid lives is denoted by d. In the context of AdS/CFT, the dual 
AdSrf + i space has d + 1 spacetime dimensions. 

Latin alphabets A,B, . . . are used to denote the d + 1 dimensional bulk indices which 
range over {r, 0,1,..., d— 1}. Lower Greek letters //, v, . . . indices range over {0, 1, ... , d— 1} 
. The co-ordinates in the bulk are denoted by X A which is often split into a radial co- 
ordinate r and x^. The dilaton field in the bulk is denoted by <5 and its value in the 
boundary is denoted by 4>. 

Our convention for the Riemann curvature tensor is fixed by the relation 

[V M ,V,]V A = -R^ x V a . (C.l) 



In Table we list the physical meaning and the definitions of various quantities used 
in the text. 



Symbol 


Definition 


Symbol 


Definition 


d 


dimensions of boundary 


Gab 
9fiu 


Bulk metric 
Boundary metric 


T 


Fluid temperature 


V 


Shear viscosity 




Energy-momentum tensor 


J s 


Entropy current 




Fluid velocity {u^u^ = — 1) 




Projection tensor, g^ v + u^u v 




Fluid acceleration, 




Fluid expansion, 




Shear strain rate, 




Fluid vorticity, 




Weyl-covariant derivative 




See Q 




Riemann tensor 


•F UV 




Gfj,u\<T 


Weyl curvature 







Table 1: Conventions used in the text, with reference to the equations defining them where appro- 
priate. 
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